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Abstract—The change in the wetted area of a packing when absorbing in water a gas having a 
high latent heat of solution has been noted in the literature. No explanation for this phenomenon 
has previously been given. Experiments have been carried out using a 1 in. diameter wetted wall 
VOL. column for absorbing ammonia from a gas stream and the water rate which will just maintain a 
continouous film on the column has been found. Below this water rate the film breaks down and 
956/57 there is a considerable reduction in the area exposed for absorption. Similar results have been 
obtained by condensing water from a hot saturated air stream on to a water film. 

The results are expressed as a difference in the surface tension, at the point where the film 
breaks, between the main bulk of the liquid and a film of liquid which has become in equilibrium 
with the main gas stream. The reduction in the surface tension of the liquid due to the absorption 
of ammonia as well as the heat of solution or condensation of the gas or vapour being absorbed 
is taken into account and also the opposite effect of the heat of vaporisation when absorbing in 
hot water. 

Ripples on the water film appear to play an important part. The use of a surface active agent 
like Teepol suppresses the ripples and prevents the breakdown of the water film until an extremely 
low water rate is reached. The breakdown of the film is also found on carbon grid packings and 
the wetted area has been measured for various flow rates. 

The surface tension of ammonia solutions at various temperatures and concentrations were 
measured and used to correlate the results obtained. 


Résumé —La littérature mentionne un changement dans la surface humide du garnissage quand 
on absorbe dans l'eau un gaz ayant une chaleur latente de dissolution élevée. On n’a pas donné 
jusqu ici d’explication de ce phénoméne. Les auteurs ont réalisé des expériences avec une colonne 
& paroi humide de 1 in. de diam., pour l'absorption de 1 NH, dans un courant gazeux, et ont ainsi 
trouvé la vitesse de l’eau qui permet d’avoir un film continu dans la colonne. Pour une vitesse 
inférieure le film se désagrége, et il en résulte une réduction considérable de la surface d’absorption. 
Des résultats analogues ont été obtenus en condensant l'eau d’un courant d’air saturé chaud sur 
film d’eau. 

Les résultats sont exprimés sous forme d'une différence dans la tension de surface, au point ot 
le film se désagrége ; différence entre la tension du corps principal du liquide et la tension du film 
de liquide qui est en équilibre avec la masse gazeuse. I] est tenu compte de la réduction dans la 
tension de surface du liquide ; réduction due 4 l'absorption de l'ammoniac aussi bien qu’a la 
chaleur de dissolution ou de condensation du gaz ou de la vapeur absorbés. 

Il est également tenu compte de l’effet opposé dd a la chaleur de vaporisation quand on absorbe 
dans l'eau chaude. 

Il semble que les rides a la surface de l'eau jouent un réle important. L’utilisation d’un agent 
tensio-actif tel que le teepol supprime les rides, évite la désagrégation du film jusqu’a de trés 
faibles débits. La brisure du film se produit également sur des garnissages du type grille de 
carbone et la surface humide a été mesurée pour différentes vitesses d’écoulement. 

Les auteurs ont mesuré la tension de surface de solutions d’ammoniac a différentes températures 
et concentrations, ce qui a permis de coordonner les résultats obtenus. 
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INTRODUCTION 
Mucnu of the design and experimental work on 
absorption is based on an assumed value of the 
wetted area of the packing. In the standard 
formula for expressing the mass transfer coefficient 
in terms of the Schmidt Number, gas and liquid 
rates, the value of the exponent of the liquid rate 
varies very considerably with the liquid rates used 
and with the gas being absorbed. Prarrt [8] gives 
an indication of the effect of absorption on the 
wetted area of a packing by giving different 
Minimum Equivalent Liquor Rates for various 
absorption systems. An example is also given of 
the effect of ammonia absorption on the water 
film on a wetted wall column which is found to 
break down prematurely. This he ascribes to the 
latent heat of solution of the ammonia but no 
further explanation is given. 

Other workers have noted that the packing of 
a tower has not been completely wetted during 
an absorption process. Moistap and Parsty [7] 
found that the wetted area of a packed tower 
varied from 10-20%, to 80-90%, when absorbing 
ethanol vapour in water at various water rates. 
They found the value of Kga depended on the 
water rate to a greater extent than was expected. 
Cuitton, Durrey and Vernon [1] found in- 
consistent results in ammonia absorption as a 
result of channelling of the liquid even though 
they had good initial liquid distribution. KowaLkeE, 
Hovucen and Warsow [5] reported the altering 
of the effective distribution of the water in a 
packed tower when absorbing ammonia from air. 
Many other workers have found inconsistent 
results when absorbing a highly soluble gas and 
have concluded that there is a change in the 
interfacial area. 

The first direct measurment of the wetted area 
of a packing was carried out by Mayo, HunTER 
and Nasu [6] using an air and water system but 
with no absorption taking place. Similar results 
were obtained by Grimiey [3] using an air and 
water system. Weisman and [11] and 
SauLmaAN and De Gourr [9] also used an air and 
water system but obtained conflicting results by 
different methods. More recently SHULMAN, 
Prou.x and ZimMERMAN [10] measured 
the mass transfer coefficient and the wetted area 
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of naphthalene Raschig rings and Beryl saddles. 
From a correlation of jp and a modified Reynolds 
number for the gas rate, the mass transfer 
coefficient (kg) for ammonia absorption was 
calculated and compared with the value of k,a 
obtained by FELiiNGerR [2]. In this way the 
effective interfacial area (a) was determined for 
ammonia absorption and it was found that it 
varied considerably with the water rate. 

Although it has been recognised for some time 
that the wetted area of a packing can be modified 
during an absorption process, little work on this 
phenomenon has been published. This investiga- 
tion was carried out to try to explain this pheno- 
menon. 


APPARATUS AND PROCEDURE 

Before the results obtained could be correlated on the 
basis of a surface tension difference, the surface tension 
of ammonia solutions in softened water was determined 
for concentrations between 0 and 13 g/l. and temperatures 
between 20 and 50°C. As expected the surface tension 
decreased with the increase in ammonia concentration. 
Jaeger’s Maximum Bubble Pressure method was used for 
the determinations. Softened water was used in the 
column to avoid the deposition of salts on the glass 
surface. Distilled water would have been preferable but 
the quantities required would be considerable. 

The initial absorption experiments were carried out in 
a lin. diameter and 4 ft long wetted wall column. To 
obtain reproducable results, great care was taken to line 
up the column vertically in order to get even distribution 
of the water on the top of the column. With no absorption 
taking place, a continuous film of water was maintained 
at water rates down to 0-05 ¢.c./sec cm. 

The gas was sampled at the inlet and outlet of the 
column by drawing a known volume of gas through N 
sulphuric acid. The solution was back titrated to determine 
the concentration of ammonia in the gas stream. A 
similar method was used for determining the concentration 
of the ammonia in the water leaving the column. The 
temperature of the gas and water entering and leaving 
the column and the humidity of the air entering the 
column were taken. Flow measurements were made by 
calibrated orifice meters. 

The hot saturated air used in one series of experiments 
was obtained by mixing steam and air before entering 
the column. Some steam was allowed to condense out 
before entering the column to ensure complete humidifi- 
cation of the air. Provision was also made for heating the 
water and air supplied to the column. 

The column was cleaned with dilute nitric acid; then 
washed with water before the start of each experiment. A 
high flow rate of softened water was then passed down 
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the column and the air and ammonia (or steam) flows 
adjusted to the required amounts. The water rate was 
reduced slowly and the condition of the film viewed 
through the glass wall. The point at which the film started 
to break at the base of the column was readily seen. If 
the water rate was reduced still further, the whole of the 
film on the column was broken and the water ran down 
one side in an irregular manner. The rest of the column 
was left dry. To reform the film, the water rate had to 
be increased considerably before a continuous film was 
formed. By careful adjustment of the water rate, the 
breakdown of the film could be controlled so that it just 
started at the base of the column. Occasionally the small 
area of column which had become dry, was rewetted but 
soon another dry patch formed at the column base. It 
was at this minimum water rate when the film had just 
started to break at the bottom of the column that the 
experiments were carried out and the results obtained. In 
all the experiments this minimum water rate was above 
that water rate which would maintain a continuous film 
on the column when no absorption was taking place. 


RESULTS 


The breakdown of the film on the column was 
initially a breaking of the surface of the film and 
consequently the surface tension between the gas 
and the liquid must play an important part. The 
first method of correlating the results was based 
on the change in the surface tension as the water 
moved down the column. From the mass transfer 
coeflicients and the surface tension results, the 
change in surface tension per unit height could 
be calculated for the base of the column where 
the film initially started to break. This, however, 
did not give any indication of the mechanism of 
the breakdown of the film. The second method, 
which is given here, was based on the difference 
in surface tension between adjacent parts of the 
film at the base of the column. Just prior to the 
film breaking, a thin film was present on the 
column wall; this film then disappeared into the 
main portion of the water as a result of its lower 
surface tension. A dry patch was then formed on 
the column and an increase in the water rate 
was required to reform the film. Now the maximum 
temperature and concentration that the thin film 
can reach is when it is in equilibrium with the gas 


stream. The difference in surface tension at the 
base of the column between the bulk of the liquid 
and the thin (or equilibrium) film was used to 
correlate these results, The surface tension of the 


The effect of absorption on the wetted area of absorption towers 


239 


liquid at the bottom of the column was found 
from the analysis of the effluent. The actual 
surface concentration and temperature will be 
slightly higher than those used but they cannot 
be readily determined. The temperature and con- 
centration of the equilibrium film was obtained 
from the inlet gas conditions and a plot similar to 
Fig. 1 was obtained for each experiment. This 
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Fic. 1. Equilibrium film conditions at base of column. 


was constructed by plotting the concentration of 
ammonia solutions at various temperatures which 
had a partial pressure of ammonia equal to the 
inlet gas concentration. Another line was obtained 
for solutions having a partial pressure of water 
vapour equal to that in the inlet gas. The inter- 
section of the two lines gave the temperature and 
ammonia concentration of the solution which was 
in equilibrium with the inlet gas to the column. 
The temperature of the equilibrium film was 
generally about 10°C and had a high ammonia 
concentration. At these concentrations and 
temperatures there is very little alteration in the 
surface tension with small temperature changes, 
and hence the surface tension values used are 
taken from the International Critical Tables. 


Concentration of ammonia varied 


The first experiments were carried out at room 
temperature using 5-50°%, ammonia in the gas 
stream and total gas velocities of 91-5 cm/sec, 
152-0 cm/sec, and 244-0 cm/sec. The water rate 
at which the film just started to break was plotted 
against the surface tension difference as shown in 
Fig. 2. In these experiments no account was taken 
of the effect of evaporation of the water into the 
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Fic. 2. Absorption of ammonia. 
A — Inlet gas velocity 91-5 cm/sec 
B - Inlet gas velocity 152-0 em/see 
C — Inlet gas velocity 244-0 cm/sec 
Inlet concentration of ammonia in air in percentage is 
shown by dotted line. 


gas stream as this only reduced the temperature 
of the water by a small amount. For a given gas 
velocity, higher water rates were required to 
maintain the film on the column as the ammonia 
gas concentration in the air was increased. 


Temperature and concentration varied 


In the second series of experiments, the tem- 
perature of the incoming water was raised 
successively from 25°C to 70°C so that there was 
considerable evaporation to counter the effect of 
the heat of solution of the ammonia. Fig. 3 shows 
the resultant heating or cooling of the water at 
the base of the column, calculated from the mass 
transfer coefficients, latent heats of solution and 
evaporation against the minimum water rate. 
Separate lines were obtained for each ammonia 
gas inlet concentration. As the inlet water 
temperature was raised, the heat lost by evapora- 
tion increased and reduced the effect of the heat 
of solution of the ammonia with a resultant drop 
in the minimum water rate. Eventually the heat 
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Iniet woter temperature increases 
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Fic. 3. Heating or cooling effect on the minimum water 
rate for various inlet ammonia concentrations in air (%). 
See Table 1 for inlet and outlet water temperatures. 


lost by evaporation exceeded the heat gained by 
absorption of the ammonia. However, the water 
film was still liable to break down due to the 
effect of the ammonia concentration on the 
surface tension. Fig. 4 gives these results plotted 


Inlet water temperature increases 
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Fic. 4. Absorption of ammonia in hot water (25-70°C). 

Inlet ammonia concentration given in percentage. Total 

gas velocity 152-0 cm/sec. See Table 1 for inlet and outlet 
water temperatures. 


on the basis of the surface tension difference 
between the main film and the equilibrium film. 
Table 1 gives the inlet and outlet temperature 
of the gas and water for these results. 
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Table 1 


Ammonia Minimum Inlet 
in Water Water 
Inlet Rate Temperature 
Gas % c.c./sec cm Cc 


| 


Outlet Heating or cooling Surface 
Water of water at base of Tension 
Temperature the column Difference 
g cal. /hr cm? ergs /cm? 


36-6 
23-9 


41-6 


170-0 
118-0 
65-6 
49-0 + 142 
55-0 44-5 3-20 


289-0 
178-0 
1-9 


92-1 


Absorption of water vapour 


To show that the concentration effect was not 
the only factor which affects the breakdown of 
the film, some experiments were made using an 
air stream heated to between 48 and 50°C and 
saturated with water vapour. The humid air 
stream was cooled by a falling film of water 
initially at 18°C and water vapour condensed 
In this 
case there was a heating of the surface but no 


and raised the temperature of the film. 


concentration effect. The film on the column was 
found to break down in a similar manner as 
before. The temperature of the equilibrium film 
was the wet bulb temperature of the gas and from 
this the surface tension of the equilibrium film 
was found. The surface tension of the main film 
was found from the temperature of the water 
leaving the column. As the difference between 
the surface tension of the main film and equili- 
brium film increased, a greater water rate was 
required to maintain the continuous film on the 
column. Fig. 5 gives the results obtained but in 
this case no effect of gas velocity, as in Fig. 2, is 
is found despite considerable variation. 


The addition of a surface active agent ( T'eepol) 
to the water supplied to the column suppressed 
the ripples and no breakdown of the film was 
observed until a very low water rate was reached 
(0-23 em/sec cm) in both the ammonia and water 
vapour absorption experiments. This low water 
rate was independent of the gas in the column. 
Furthermore, no breakdown of the film was 
observed when evaporating water from a falling 
film. 
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Surface tension difference at 
moment of breakdown of the film ——»> ergs/cm* 


Fig. 5. 


Absorption of water vapour in cold water. 
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5-35 0-34 39-5 37-3 - 23-9 475 
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9-5 0-439 | 43-0 41-7 5-50 
10-5 0-306 49-8 46-9 - 3-87 
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20-0 1-110 34-4 
19-7 1-060 0-5 
19-6 0-891 419-0 
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27-6 1-450 27-6 32-0 + 11-28 
28-0 1-450 0-9 43-8 7-78 
28-5 1-350 54-8 56-3 5e25 
27-5 0-778 68-6 65-5 2-00 
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Absorption on serrated carbon grids 

Some experiments were also carried out to 
determine the wetted area of an industrial type 
of packing for various gas concentrations using 
lin. carbon grid, }in. thick and 3gin. wide, 
marked out in squares. 

An even distribution of water was maintained 
to the grid and at various gas and water rates the 
area wetted by the water was noted. With no 
absorption the grid was completely wetted at a 
water rate of 0-2 c.c./sec. The results are shown 
in Fig. 6. No change in wetted area was found 
with varying velocities but this may have been 
due to the low velocities used (13-2 cm/sec 
maximum). In these results it was particularly 
noticeable that the water flowed over the grid 
in a parabolic manner. 


2 100 


Wetted crea 

| 


a 


° 1 2 3 5 
Water rote cc /sec 
Fic. 6. Wetting of a carbon grid while absorbing ammonia 


or water vapour. Ammonia concentration in air in %. 
Inlet saturated air temperature in °C. 


The water fell from the drip point of the grid 
on to the top of the grid below and the majority 
of the water spilt over on to its face. The main 
film thus started from a point source and spread 
out giving a parabolic pattern. The remaining 
water ran along the top of the grid and overflowed 
on to the face. This secondary water formed a 


small thin film between the parabolic films beneath 
the drip points. The parabolic film was compara- 
tively thick compared with the thin film between 
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but it was the latter film which broke first and 
formed a dry patch when absorption was taking 
place. In some cases, particularly at the higher 
ammonia gas concentrations, the main parabolic 
film also broke down when absorption was taking 
place. 


Discussion OF RESULTS 


The presence of ripples appear to play an import- 
ant part in the breaking down of the film. When 
the ripples are suppressed by a surface active 
agent, the film does not break down while 
absorbing ammonia or water vapour until a very 
low water rate is reached. This low water rate 
is independent of the absorption of the ammonia 
or water vapour. 

Grmey [8] found that ripples start at a 
Reynolds Number of 25, reach a maximum of 
intensity at 800 and cease at a Reynolds Number 
of about 1,000, after which turbulence appears 
to set in. Plotting the results using the Reynolds 
Number for the minimum water rate, Fig. 7 is 
obtained. The gradient of curve B decreases as 
the water rate approaches a Reynolds Number 
of 1,000. Curves C and D follow a similar line 
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A — Absorption of water vapour 
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but it was not possible to reach the higher Rey- 
nolds Number in the experiments. Curve E 
reaches a maximum at a Reynolds Number of 
1,080. At the higher water rates it was often 
found that the dry patch formed at the base of 
the column was rewetted, only to be followed by 
another dry patch. At high water rates the 
turbulence in the water film appeared to reduce 
the effect of the ripples and rewet the column 
whenever a dry patch was formed. 

The presence of ripples on the column gives 
rise to an unevenness in the film. In the case of 
ammonia absorption, the ammonia diffuses 
through the air to the surface of the water and 
dissolves but the diffusion in the water layer is 
greater where there is a thick film at the crest 
of a ripple and hence larger volume of water, 
than at the trough where the film is thin. The 
trough will then have a higher concentration of 
ammonia than the crest and the maximum 
concentration it can reach is when it is in equili- 
brium with the gas stream. The temperature of 
the trough will also be determined by the equili- 
brium conditions. As a result of the equilibrium 
conditions for the trough, the surface tension will 
be lower than that of the crest or thick film. This 
difference is sufficient to draw the thin film into 
the thick film leaving a dry patch on the bottom 
of the column. An increase in the water rate 
is then required to rewet the dry part of the 
column. 

As the temperature of the inlet water is 
increased during ammonia absorption, the 
evaporation of the water begins to play an 
important role and is a counter-effect to the heat 
of solution. The equilibrium film remains at 
substantially the same temperature and con- 
centration for a given inlet gas concentration as 
the inlet water temperature is increased. The 
temperature of the main bulk of the liquid 
at the bottom of the column will depend on the 
amount of absorption and evaporation that is 
taking place. At low inlet water temperatures 
there is a resultant heating of the liquid but, 
as the inlet water temperature is increased, more 
water is vaporised and less ammonia is absorbed. 
At the highest inlet water temperature used the 
heat of vaporisation is greater than the heat of 


solution of the ammonia and the outlet tempera- 
ture of the liquid is reduced. The surface tension 
of the main bulk of the liquid at the base of the 
column is affected by the temperature as well as 
the ammonia concentration, the net result being 
a decrease as the inlet water temperature is 
increased. Since the surface tension of the 
equilibrium film is substantially constant then 
the surface tension difference at the base of the 
column is decreased and a lower minimum 
wetting rate is required for the column. Even 
though the heat lost by evaporation can exceed 
the heat gained by the heat of solution of ammonia 
the film on the column can still break down owing 
to the concentration effect, the surface tension 
of the bulk of the liquid being slightly higher 
than the surface tension of the equilibrium film. 

The low water rate at which Kanryka and 
Hinckuierr [4] were able to operate their 
adiabatic hydrogen chloride absorber was probably 
due to the large amount of evaporation of water 
that was taking place. The heat lost by evaporation 
counteracted the heat gained by the solution of 
hydrogen chloride with the effect that the surface 
tension difference was reduced. A much lower 
water rate than normally used in hydrogen chloride 
absorption was possible while still maintaining 
efficient wetting of the packing. 

In the case of the condensation of water vapour 
on the falling film, there was no concentration 
effect. However, the thin film in equilibrium 
with the gas stream would have the same tem- 
perature as the gas stream. The main part of the 
film would have a lower temperature than the 
humid gas stream and hence a higher surface 
tension than the equilibrium film. The surface 
tension difference was sufficient to form a dry 
patch on the column. 

When there is evaporation from the water film 
into an inert gas stream, no premature breakdown 
of the film is observed. 

In the case of carbon grids, the break comes 
first between the patabolic water films below the 
drip points but at high ammonia gas concentra- 
tions the main parabolic films often break down. 
No effect of the velocity of the gas was found in 
this case but it was probably due to the very 
low values obtainable. Alternatively, it may be 


248 


. 
) 
Boss 
y 
_ 
qa 


J. Bonp and M. B. Donatp 


due to the different flow patterns obtained on 
the grids where thick and thin films are much 
more pronounced compared with wetted wall 
columns. 


CONCLUSIONS 


The absorption of ammonia or water vapour 
on a continuous film of water with a rippled 
surface can result in its premature breakdown 
reduce the efliciency of the 


and hence 


column. 
The breakdown can be effected by concen- 
tration as well as temperature effects. 


(3) An increase in the water temperature will 
increase the heat lost by evaporation. This 
will counteract the heat gained by absorption 
and hence reduce the water rate at which the 
film breaks. 

(4) The breakdown of the film can be prevented 


by using a surface active agent. 
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The quaternary system: benzene—carbon tetrachloride—acetic acid— 
water at 25°C 
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Abstract—Apparatus is described for the accurate determination of liquid — liquid equilibrium 
in ternary and quaternary systems by the synthetic method; and in particular, a microburette 
for use in cloud point titrations, and a flask which enables the weights of the equilibrium phases 
to be determined directly. Tie lines may then be located by the “ lever rule,” more accurately 
and conveniently than in the method of Orumer, Wuire and Trurcer [14]. Results are 
presented for the system benzene —- carbon tetrachloride — acetic acid —- water at 25°C, and 
relationships between the equilibrium concentrations in the quaternary and the ternary systems 


VOL, are discussed. 


956/57 Résumé—Les auteurs décrivent un appareillage pour la détermination précise de léquilibre 
liquide-liquide dans les systémes ternaires et quaternaires par la méthode synthétique. Ils 
décrivent en particulier une microburette utilisable dans les titrages en point de trouble et un flacon 
qui permet de déterminer directement les poids des phases en équilibre. Les lignes de raccorde- 
ment peuvent étre alors déterminées par l'utilisation de la * courbe de distribution * plus pratique 
et plus adéqu ate que dans la méthode d’Oramer, Warre, et Truecer. Les auteurs donnent des 
résultats pour le systeme benzéne — tétrachlorure de carbone — acide acétique eau a 25°C et 
discutent les relations entre les concentrations d'équilibre dans les systemes térnaires et 


quaternaires, 


INTRODUCTION system illustrated in Fig. 1(a), in which one 
pair of components is partially miscible. 
Cuanc and Mouton in their study of several 
quaternary systems of the type illustrated in 
Fig. 1(b), where one component is partially 


EquiLtisrium relationships in four component 
two phase liquid systems are required in the 
design of solvent extraction processes involving 
complex solvents. As it is not often possible 
to calculate such equilibria theoretically, and 
practical determinations are difficult and tedious, 
a first step is to relate the equilibrium concentra- 
tions in the quaternary systems to those in the 
component ternaries, about which more is 
generally known, or which may be determined 
comparatively easily. Quaternary systems have 
been studied for this purpose by BRraNcKER, 
Hunter and Nasu [4] and, more recently, by 
CuancG and Mouton [7]. The former postulated 
a comprehensive set of relations between the 
quaternary and ternary equilibrium concentra- 
tions, which will be discussed below. These 
relations were applicable to the type of quaternary Fic. 1. 


miscible with two others, found a relationship 


*Present address: The Distillers Co. Ltd., Research and Development Department, Great Burgh, Epsom, Surrey. 
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between the quaternary equilibrium concentra- 
tions and one of the ternary systems, viz. the 
double non-consolute pair ternary BCD of Fig. 
1(b). It was decided to investigate a number of 
quaternary systems in order to examine further 
the existence of such relations. The apparatus 
and experimental techniques evolved for this 
purpose, and the results obtained for a quaternary 
system of the type shown in Fig. 1 (b), namely 
benzene — carbon tetrachloride — acetic acid — 
water, are presented here. 


EXPERIMENTAL METHODS 


The apparatus and techniques employed will be 
described briefly, as they are novel in some respects, and 
appear to have advantages over the more traditional 
methods. It was necessary that methods be selected 
which could be used with the widest possible range of 
systems, and that they be sufficiently accurate to discri- 
minate between competing correlations (this point has 
been discussed in detail by Hann [9)). 

The first of these criteria eliminated direct analysis 
of the equilibrium phases, and the alternative “* synthetic ™ 
method was adopted. This two step method requires 
(i) the determination of the saturation isotherm by 
titration, and (ii) the location of tie lines, using information 
obtained in (i) to identify the tie line compositions. The 
usual method of tie line location involves in an mn component 
system analysis of the phases for n-2 components, or 
measurement of n-2 physical properties. 

A more elegant method of tie line location is based on 
the “ lever rule,” which states that the point representing 
the initial mixture lies on a line joining the points represent- 
ing the compositions of the two equilibrium phases, and 
divides this line in the inverse ratio of the weights of the 
two phases. If the initial mixture composition is known, 
measurement of this ratio is then sufficient to determine 
the tie line in a ternary system provided the saturation 
curve is available. This method was applied in early 
work [3, 11], but in both cases invalidating assumptions 
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were made. Valid determinations of ternary tie lines 
were made by Orumer, Wuaire and Truecer [14]; 
they made up mixtures in a graduated cylinder, so 
that, after attaining equilibrium and allowing the 
phases to settle, phase volumes could be read off. 
Density measurements were required to convert the 
volumes to weights. A graphical method for locating tie 
lines from such data was given later [13]. 

This method when applied to a quaternary system 
requires one additional measurement, which can con- 
veniently be a physical property of one of the saturated 
phases. In Fig. 2 the saturation frustrum of the quaternary 
system type illustrated in Fig. 1 (b) has been reproduced. 
DE is a quaternary tie line which passes through the 
overall or initial mixture composition point M; and 
points D and E lying on the saturation surface of the 
frustrum represent compositions of quaternary phases 
in equilibrium. A convenient physical property of the 
phase E has been ascertained to be of value z. The weight 
of the two phases D and E has been determined, while 
the overall mixture composition M is also known. The 
point E must lie on the line FG of constant physical 
property (value z), which line lies on the saturation surface 
and has been obtained in the course of determining that 
surface. Further, the quaternary tie line DE must be 
so located in space that the ratio DM/ME equals the 
ratio (weight of phase E)/(weight of phase D), and it 
may then be found by trial and error. This procedure 
was used in the present work. The method differed in 
practice from OruMer’s in that the phases were separated 
and weighed directly, rather than that the weights were 
inferred from phase volume and density determinations. 


APPARATUS 


In the cloud point titrations (by which the 
saturation isotherm was established) and the 
tie line determinations agitation of the solutions 
was provided by a shaker mechanism which 
rocked the solution flask through an are of 10 
at speeds which were continuously variable from 
90 to 900 oscillations per minute. 

The number of points to be determined in the 
cloud point titrations limited the quantities for 
each point to such that final volumes of about 
20 ml would be obtained. It was then desirable 
to add titrant in increments of less than 0-002 ml, 
and hence a microburette was required. Although 
a large number of designs of such instruments 
have been published (reviewed by Kirk [10)), 
none appeared to be completely satisfactory 
for the present purpose. After extensive trials 
of several types, the burette shown in Fig. 3 was 
constructed. 
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The quaternary system : 


This consisted of a stainless steel reservoir, 
a glass vessel containing the titrant liquid, and a 
flexible connecting unit, to convey the titrant 
from the stationary burette to the shaking flask. 
The dimensions of the reservoir gave a convenient 
displacement, 8-5 ml. A rubber “ piston ring,” 
as illustrated, was found necessary to prevent 
mercury leakage. 24 BG stainless steel tubing 
was used for the connecting unit, and fused 
directly into Imm Pyrex capillary 
In operation the burette was filled with sufficient 
clean mercury to cover the glass-metal joint, 


tubing. 


ensuring that the titrant was kept to the titrant 
vessel and the connecting unit, and hence only 
these needed thorough cleaning when changing 
titrant. After filling the burette with titrant the 
tip of the steel tube was arranged to be within 
a half inch of the solution surface, so that during 
shaking titrant was washed off the tube tip, and 
hence could be added in increments well below 
the desired figure. 

The burette gave easy and positive control 
of the rate of titration, and could be mounted 
conveniently on the thermostat tank. The 
titrant was protected from atmospheric dirt 
and moisture, and from evaporation, and had no 
tendency to leak into the solution to be titrated. 
Construction in three separate parts made for 
easy cleaning and filling. 

The solution to be titrated was contained in 


ff Connecting unit 
Imm capillary 


24 BG SS tube Bin. long 


A7 cone 


B24 quickfit 


50 mi. titrating flask 


Fic. 3. Titration assembly. 
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benzene carbon tetrachloride acetic acid-— water at 25°C 


the special stoppered flask, also illustrated in 
Fig. 3, in which a volume of 20-25 ml of liquid 
could be safely shaken at speeds of up to 900 
oscillations per minute without loss of liquid. 
For filling and cleaning the stopper was removed, 
but during titration only the cap on the stopper 
capillary tube was removed and the steel tube of 
the microburette inserted through the capillary 
tubing. Even during a lengthy titration there was 
little danger of incurring significant losses due to 
evaporation. 

The apparatus used for the location of tie 
lines is shown in Fig. 4. The bulb was intended 
to contain about 40 ml of liquid, and had a total 
volume of about 100 ml. It was filled through the 
bent neck, which prevented any loss of liquid 
The 
diameter of the capillary tube was found to be 


due to splashing onto ground surfaces. 


critical, too small a diameter causing excessive 
rise of liquid, and hence loss, while too large a 
diameter made phase separation difficult or 
impossible, due to the formation of bubbles of 
one phase in the tube. The optimum diameter 
yas found to depend on the surface tensions 
of the solutions, and that shown (1 mm) was 
satisfactory the 
but more recent work with this apparatus on a 


found for system studied, 


non-aqueous system required smaller diameter 
tubing. The bore of the capillary should also be 


A7 cone BW cone 


Imm capillary 


4 in. tubing 


Fic. 4. 


Equilibrium flask. 
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uniform so that small quantities of a phase may 
be estimated. Both neck and capillary tubes 
were normally closed by mating ground caps. 

In making a particular determination, the 
requisite quantities of components were weighed 
into the vessel, which was then placed in the 
thermostat and shaken at about 700 oscillations 
per minute. Shaking was stopped at frequent 
intervals, the phases allowed to settle out, and 
liquid displaced from the capillary tube by 
applying pressure from a dry air supply, so as 
to bring this liquid to equilibrium with the bulk 
of the solution. As this process produced marked 
turbidity until equilibrium was established, this 
event could be confirmed visually. After an hour 
the shaker speed was reduced to 150 oscillations 
per minute, producing gentle agitation. This 
enabled the phases to settle, while changing the 
interface, and helped to remove bubbles of the 
lower phase remaining in the upper surface of 
the upper phase. Any such bubbles still remain- 
ing could usually be removed by gentle twirling 
of the vessel after the bulk of the lower phase 
has been separated. Shaking at the lower speed 
continued for at least two hours, and was followed 
by settling for at least one hour. The lower phase 
was then sampled, if required, by pushing liquid 
out through the capillary tube by applying air 
pressure to the neck, followed by removal of the 
bulk of the lower phase in the same way. The 
last portions of the lower phase were removed 
by holding a folded piece of filter paper over the 
end of the capillary, and in this way a neat 
separation of the phases could be attained. 
The vessel was weighed after phase separation, 
giving the weights of the upper and lower phases. 

Densities were determined using 10 ml Sprengel 
—Ostwald type pycnometers. Samples were 
pushed into them from the titration or equilibrium 
flasks through special sampling tubes, which 
prevented the samples from coming into contact 
with the outside atmosphere. 

All weights were determined to within 

0-0001 g. Flasks and pycnometers were 
weighed against tares, and in all cases the usual 
precautions were taken to ensure that glassware 
which was wetted during an experiment would 
have a reproducible film of moisture on it. 
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The cloud point titrations and phase separations 
were carried out in a water thermostat controlled 
to 25-00 + 0-02°C at all points. Three sides of 
the thermostat were made of plate glass to permit 
illumination and observation of the solutions. 


MATERIALS 


Water was redistilled from laboratory distilled 
water, to which alkaline potassium permanganate 
had been added, in an all Pyrex apparatus of the 
type standard for conductivity water, using an 
entrainment trap [15], the middle third of any 
batch being collected. The product was stored 
protected from atmospheric carbon dioxide by 
self indicating soda asbestos. 

Properties: Conductivity: 0-7! 10° mho/em 
at 25.00°C, 

Acetic Acid, A.R. quality, was dried by partial 
freezing until the freezing point exceeded 16°C, 
A.R. acetic anhydride was added, 20°, in excess 
of that required to convert the water content 
(calculated from the depression of the freezing 
point) to acetic acid, and the acid allowed to 
stand for some days before fractionation in a 
column packed with 3 mm Fenske helices. After 
refluxing for 15 hours, 20% of the charge was 
rejected at a reflux ratio of better than 10: 1, 
and 40°, collected at a ratio of better than 
4:1. Norise in boiling point was observed during 
the collection. Acid was prepared on two 
occasions, a 2ft column being used for batch 
1 and a 4ft column for batch 2. 


Properties : Batch 1 Batch 2 
B.P. (760 mm.) 117-9 17-9 C 
Density, 25°C 1.0438 1-0438¢ /m] 
F.P. 16-52 16-55°C 
Water content 0-043 0-028°%, 


Benzene was also prepared in two batches, 
from thiophene free material. For batch 1 
benzene was refluxed and fractionated over 
sodium in the shorter column, twice recrystallized, 
and again refluxed and fractionated over sodium, 
using the same procedure as for the acid. 
For batch 2 the benzene was washed with 
chlorsulphonic acid, neutralized and washed with 
water, refluxed and fractionated over sodium 
in the 4ft column, and recrystallized twice. 


oe 
VOl 
6 
195¢ 
4 
A 
d 
|| 


VOL, 
6 
956/57 


The quaternary system : 


Properties : Batch 1 Batch 2 
B.P. (760 mm) 80-1 80.1°C 
Density, 25°C 0.8732 0.8736 g/ml 
F.P. 5-40 5-47°C 


Carbon tetrachloride, technical grade, was 
washed in half gallon lots with 100-200 ml 
consecutive of alcoholic potassium 
hydroxide, distilled water, sulphuric acid, water, 
alkaline potassium permanganate, and _ finally 
water. The product was allowed to stand over 
calcium chloride for several days, and was then 
refluxed and fractionated over dried anhydrous 
A.R. sodium carbonate, as for the acid. 

Properties : 

B.P. (760 mm) 
Density, 25°C 


washes 


76-7°C 

1-5846 g/ml 

The acetic acid, benzene and carbon tetrach- 
loride were stored in storage bottles fitted with 
delivery burettes, both bottles and burettes 
being kept under an atmosphere of nitrogen 
which had been dried over anhydrous magnesium 
perchlorate, and from which carbon dioxide 
had been 
asbestos. 


removed by  self-indicating soda 


THE QUATERNARY ISOTHERM 


The quaternary saturation surface was explored 
by taking sections through the acetic acid and 
water apices, and finding the saturation curve 
for each section. Each section could be regarded 
as a ternary system for this purpose, of com- 
ponents acid, water, and a_ benzene — carbon 
tetrachloride mixture. For each section a 
sufficient quantity of the benzene — carbon 
tetrachloride mixture was made up by weight, 
and stored under conditions making a significant 
change in composition unlikely. The quaternary 
isotherm was determined as nine such sections, 
five of these being complete sections of 16 experi- 
mental points each, and four sections of eight 
points only on the water lean side of the isotherm, 
making 112 experimental points in all. All 
points were found by cloud point titration, using 
the apparatus described. For all points in the 
8-point sections, and for 12 of the 16 points in 
each of the complete sections, water was the 
titrant, as preliminary experiments had shown 
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benzene — carbon tetrachloride — acetic acid- water at 25°C 


that where both could be used, titration with 
water gave better cloud points than titration 
with the benezene — carbon tetrachloride mixture. 

The results are shown in Table 1. In calculating 


the observed compositions, the water content 
of the acetic acid and the benzene, found from 
their freezing points, was taken into account. 
The specific volumes found for the saturated 
solutions are also shown. Two series of results 
are given for the benzene — acetic acid — water 
ternary, the first being carried out with benzene 
batch 1, and repeat determinations being made 
with the purer batch 2, giving the same solubili- 
ties as before, within experimental error, but 
slightly differing specific volumes. 

As the solubilities in the various sections are 
similar, individual results were corrected for 
over-and under-titration by comparison with 
the other results, and three charts were made up 
for the purpose. In the first, the results for the 
nonaqueous phase with acetic acid concentrations 
up to 30% were plotted as weight per cent 
water against weight per cent acid on a scale of 
10 cm to 1% water and 1 cm to 1% acid. Results 
for the aqueous phase were similarly plotted, 
as benzene + carbon tetrachloride (10 cm to 1°) 
against acetic acid (0-5 cm to 1%), covering the 
range to 54% acid. The large changes in composi- 
tion of all components over the remaining, 
central range made a similar arithmetic plot 
unsuitable, and Bancrort’s expression for the 
isotherm [1] was used, log (water/acetic acid) 
being plotted against log (benzene + carbon 
tetrachloride) /acetic acid on a 60 cm 100 cm 
chart. The necessary binary data were taken 
from SEIDELL [20]. 

For most of the points over-or under-titration 
had been recorded by noting whether a titrated 
solution became cloudy after heating to clearness 
and bringing back to 25°C. Using this information 
smooth curves were drawn on the above charts 
to represent the saturation curves of the sections. 
For points not lying on their appropriate curves, 
the composition of the “ adjusted” saturated 
solution was calcu'ated from the curve. The 
“ adjusted” values so obtained are given in 
Table 1. It will be noted that while most of the 


corrections amount to changes in composition 
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Table 1. Cloud point titration results for the quaternary isotherm Benzene (B)— Carbon Tetrachloride 
(C) — Acetic Acid (A) — Water (W) at 25°C 


| Observed Adjusted 
| Specific | 
Point) Bic Bac A | Volume B | Cc A 
No. Wt. % we % | | % | 25) | | We % | | We % 
ml/g at 4 
1 0 4-82 | 4-96 o22 | 11ss1s | 94-86 94-86 | 4-97 0-17 
2 Batch (1) 89-52 | 1013 | O85 | 1-13057 89-54 89-54 10-13 0-33 
3 Benzene 82-40 16-93 | 0-67 111878 | 82-43 82-43 16-93 | 0-64 
74-41 | 24-39 1-19 101502 | 7439 | 7439 | | 2438 | 1-22 
5 | 67-83 30-29 189 | 109261 | 67-84 67-84 | 30-29 | 1-88 
6 59-63 | 37-26 3-11 1-07612 59-61 | 5961 | 37-25 S14 
7 | | 49-60 45-16 524 | 105513 | 49-60 19-60 | $5-16 5-24 
. | 35-87 | 54-98 9-16 1-02342 35-95 | 35°95 | 55-10 8-06 
9 26-56 61-05 12-39 | 100263 | 26-57 26-57 6108 | 1235 
| 16-16 | 6572 1811 | | 1613 | 1613 65:59 | 18-27 
| 11-61 65-98 2242 | 0-96699 11-38 | 1138 65-83 22-60 
12 +71 59-51 35-78 0-95359 +-70 59-45 35-85 
13 | 2-82 53-54 3-64 | O-95232 | 2-84 2-84 | 5B-5s $3-63 
4 | 1-57 $4-56 53-87 | O-95554 1-57 1-57 $4-56 53-87 
5 | 0-67 29-70 69-6: 0-96703 0-67 0-67 29-70 =| 69-63 
i 0-98325 0 | O40 14-89 84-72 


17 0 4-76 oo 0-34 1-13773 94-02 94-92 4-91 0-17 
18 | Batch (2) 8o-62 10-03 0-36 112974 | 89-64 89-64 10-03 0-33 
19 Benzene | 82-24 17-11 066 | | 82-24 82-24 711 (O65 
20 75-08 23-81 1-12 110544 | 75-04 75-04 | 2380 | 1-16 
21 68-08 30-09 183 09251 68-07 6807 | 30-08 1-85 
22 5006 | 4681 | 1-05556 5006 =| 5006 | 5-31 
23 36-78 ee 89 1-02550 36-84 36-84 54-46 8-67 
24 12-412 | 0-22 1-07723 94-72 82-96 11-76 5-11 0-17 
25 80-07 | 10-01 0-31 107268 80-67 78-54 11-13 | 
26 82-18 17-18 0-64 1-06542 82-18 71-98 | 10-20 | 17-18 0-64 
27 74-79 | 24-03 118 | 105740 74-80 65-51 | 928 | 2408 1-17 
28 67-40 «| 30-67 1-93 1-04817 67-40 59-03 8-37 30-67 | 1-93 
29 59-87 87-05 3-08 | 103771 59-87 5244 | 743 3705 3-08 
30 49-06 44-85 5-20 1-02331 50-00 43-79 6-21 44-58 5-13 
31 | 36-58 54-51 | 8-91 1-00182 | 36-66 32-11 | 4°55 | 54-64 8-70 
32 | 24-953 04-71 5-06 0-22 101610 | 94-77 71-12 | 23-65 | 5-06 0-16 
33 0-34 | 1-01508 809-90 | 67-47 22-43 9-80 0-29 
: S219 | 17-16 | 0-66 101250 | 8221 | 61-70 | 2051 | 
35 7488 | 2305 |) 1-17 100804 | 7490 | 56-21 | 18-69 23-05 | 115 
36 | 6728 | 30-74 | 1-98 100429 | 67-31 | 50-51 16-80 | 30-76 1-93 
37 | 60-34 3664 | 3-02 090033 | 6035 | 4520 | 15-06 36-65 3-00 
sa | | 49-92 482 5-27 0-99057 50-03 | 37-55 | 12-48 44-91 5-07 
39 | 26-37 | 54-07 8-06 0-97766 | 36-48 27-38 9-10 54-83 8-69 
26-53 61-15 12-32 06706 26-56 19-03 6-63 61-23 12-21 
4! 16-57 65-65 i7-78 | O@5714 16-57 12-44 $15 65-65 17-78 
‘ 24-953 ; 1189 | 65-09 22-12 0-95214 | 1189 | 892 207 | 65-99 22-12 
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benzene — carbon tetrachloride - acetic acid —- water at 25°C 


The quaternary system : 


Table 1 (continued) 


Specific | 
: Point} B+C B+cCc A w Volume | B+C B | ¢ A w 
No. we % Wt % | Wt % | WE % | as | Wt % | Wt % | Wt % wy | we 
‘ | | 
i 43 | 4-82 59-71 35-47 0-94735 4-83 3-62 1-21 59-90 35-27 
: 44 2-87 53-46 43-67 0-94885 2-75 2-06 0-69 53-53 43-72 
45 | 1-58 44-61 53-81 095392 1-53 1-15 0-38 44-63 53-84 
46 0-76 29-84 69-41 096637 0-65 0-49 0-16 29-87 69-49 
47 0-46 14-69 84-84 0-98335 037 | 0-28 0-09 14-70 84-92 
48 87-436 94-79 5-01 0-19 095542 94-83 59-33 35-50 | 5-01 0-15 
49 89-84 9-88 0-28 0-95775 89-84 56-21 33-63 9-88 0-28 
50 82-45 16-97 0-58 095953 82-44 51-58 30-68 16-97 0-59 
VOL. 51 74-91 23-96 1-13 096073 74-90 46-86 28-04 23-96 1-14 
6 52 67-44 30-60 1-97 096094 67-48 42-22 | 25-26 30-62 1-91 
956/57 53 59-84 37-00 3-15 0-96010 59-90 37-48 22-42 37-03 3-06 
47-39 46-68 5-93 . | 47-58 29-78 | 17-80 46-86 5-56 
55 36-34 54-77 8-88 0-95401 | 36-43 22-79 13-63 54-90 8-66 
56 49-965 94-67 5-10 0-22 089426 94-75 47-41 7-34 5-10 0-14 
57 89-93 9-80 0-27 0-89928 89-93 45-00 44-93 9-80 0-27 
58 82-65 16-7 0-57 0-90618 82-66 41-36 41-30 16-78 0-56 
59 | 75-12 23-82 1-07 0-91209 75-09 37-57 37-52 23-81 1-11 
60 67-76 30-31 1-93 091684 | 67-80 33-92 | 33-88 30-33 1-87 
61 | | 60°70 36-29 3-01 092065 60-76 30-40 | 30-36 36-35 2-91 
62 50-18 44-62 5-25 092530 50-26 25-15 25-11 44-74 5-00 
63 | 36-83 54-54 8-62 0-93004 36-89 18-46 18-43 54-63 8-48 
64 26-85 61-19 11-96 0-93316 26-85 13-43 13-42 61-19 11-96 
65 | 16-95 65-69 17-37 0-93534 16-95 848 | 847 65-69 17-37 
66 | 1152 66-08 22-40 | 0-93678 11-52 5-76 5-76 66-08 22-40 
67 §-22 60-7: 34-05  0-94088 5-23 2-62 2-61 60-81 33-96 
68 2-79 53-69 43-52 094547 2-79 140 | 1-39 53-69 43-52 
69 1-52 44-30 54-18 095247 1-46 073 | O73 44-33 54-21 
70 0-76 29-68 69-56 096599 0-62 0-31 0-31 29-72 69-66 
71 0-44 14-91 84-65 098288 0-35 0-18 0-17 14-92 84-73 
72 62-460 94-93 4-96 O-11 083354 94-92 35-63 59-29 4-96 0-12 
7% 89-93 9-92 0-16 084169 89-85 33-73 56-12 9-91 0-25 
74 82-45 17-00 0-54 085320 82-44 30-95 51-49 1700 | 0-55 
75 74-67 24-20 1-13 086433 74-68 28-03 46-65 24-20 1-12 
76 67-68 30-39 1-93 087319 67-71 25-42 42-29 30-41 1-88 
17 59-85 36-97 3-18 0-88215 59-04 22-50 | 87-44 37-03 3-03 
78 49-78 | 44-92 5-30 089285 49-93 18-74 | 13-19 45-07 5-01 
79 | 86-73 | 54-75 8-52 090612 36-79 | 13-81 | 22-98 le 
80 75-216 04-87 502 | 0-10 0-77094 94-87 | 23-51 | 71-36 502 | O10 
81 89-98 9-83 0-19 O-78244 89-95 22-29 | 67-66 983 | O22 
82 82-51 16-98 | 0-51 0-79889 82-50 20-45 62-05 16-98 | 0-52 
3 75216 | 74-90 2402 | 1-08 081433 74:91 | 18-57 | 56-34 2402 | 107 
84 | 67-68 30:39 | 192 0-82837 6771 | 16-78 | 50-93 30-41 | 1-87 
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Table 1 (continued) 


100-000 


60-05 
50-07 
36°85 
26-91 

700 
11-92 
4-96 
2-83 
142 
0-63 
0-56 


95-00 
89-92 
83-07 
74°97 
67-60 
59-98 
50-23 
36-96 


95-07 
91-10 
84-76 
80-20 
75-40 
68-91 


Observed 


36-83 
44-80 
54°89 
61-55 
66-08 
66-49 
60-66 
54-36 
44 
29-62 
14-80 


404 

9-92 
16-49 
23-08 
30-51 
36-95 
44-82 
55-08 


8-78 
14-90 
19-20 
23-69 
20-41 
43°26 
55°48 
62-49 
66-86 
67-17 


61-78 


53-65 


44-51 


20-85 
14-97 


3-12 
5-13 
8-27 
11-53 
16-91 
21-59 
34-38 
2-81 
54-24 
69-75 


84-64 


O-O6 
O16 
Ot 

105 

3-07 


7-96 


0-07 
0-12 
0-60 
0-91 
1-68 
4-35 
763 
1144 
16-78 


21-10 


| S343 
43-76 
54-32 


69-57 
84-71 


| 


0-84197 
0-85838 
O-88142 
O-89751 
0-91282 
0-92094 
0-93485 
0-94157 
0-95096 
0-96546 
0- 98284 


0-70938 
O-72440 
0-74372 
0-76546 
0-78441 
0-80264 
0-82699 
0- 85686 
0-64915 
0-66349 
0-68553 
0-70085 
0-71652 
O-7T3701 
0-78733 
O-835635 
0-86422 


| 89272 
| 


90612 


| 


0-939044 
0-95006 
0-96408 
0-98263 


| 


60-10 
50-20 
36-88 
26-91 


14-90 

12-44 
9-14 
6-67 
$21 
2-95 
1-23 
0-71 
0-35 
O15 
0-08 


11-61 
10-99 
10-15 
9-17 
8°26 
7-34 
6-15 
4-52 


Adjusted 


45°20 | 36-86 
37-76 44-92 
27-74 54-93 
20-24 61°55 
12-79 66-08 
8-97 66-49 
3-73 60-66 
2-16 54-34 
1-07 
0-45 29-63 
14-83 


‘ 404 
78-89 0-92 
72-00 16-49 
65-81 23-98 
50-35 | 30- 
52-70 | 
+18 | 446-91 
32-45 | 55-10 


95-07 4-86 
91-07 8-78 
84-75 14-00 
80-20 19-20 
75-35 23-68 
68-90 20-4 
52°45 43°28 
46-02 55°39 


521 61-87 
2-5 53-62 
1-39 
0-58 20-85 
0-30 14-07 


| 34-38 
| 42-79 


54-24 
69-77 
84-85 


0-09 
| O20 

0-46 
1-04 
1-88 
2-98 
4-77 


7-93 


of less than 0-1°%, a few are higher and lie largely 
0 get) 


in one region, near the plait points. 


TERNARY AND 


methods described above, 
acid — water 


benzene — acetic 


QUATERNARY 


ternary 


TIE 


LINES 


Tables 2 and 3 show the data, obtained by the 
for six tie lines in the 
system, 


nary 
shown in 
titration, generally in duplicate, of samples with 
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nine lines 
acid — water ternary, and 23 lines in the quater- 
Acetic 
some cases ; 


system. 


in 


the 


carbon 


acid 


they were obtained by 


standardized N/5 barium hydroxide, using bromo- 
thymol blue as indicator [4]. 


tetrachloride 


concentrations 


acetic 


are 


| | | 
: Specific 
Point) B+C | B+C 4 Ww |B+C| B 7 
No. wey | wey | | % | wey | | | | we. % ‘ 
85 i | | | | | | 3-04 : 
86 4-88 
88 | | 11-53 
90 | | 92 21-59 
91 | | 96 
92 2-87 
a 93 | | | 1-42 
95 | | vol 
| | | 195€ 
89-88 
4 98 | 83-05 | | 
| | 74-98 | | 
100 | | 6761 | 
101 | | | 60-04 | 
4 102 | 50-33 | ic 
4 103 | | 36-97 | | 
104 | | 95-07 | 007 
105 | 91-07 | 
106 | | 84-75 0-35 
107 | 80-20 0-60 2 
4 108 | | 75-35 | | 0-97 
109 68-00 | 
4 110 | 5241 | 52-43 | 4.99 eg 
| 36-95 36-92 
| 26-47 | 26-45 | 26-45 | 6245 | 1110 
118 16-36 16-39 1639 | 6TOL | 16-60 
114 11-73 — 11-73 1-73 | 67-17 21-10 
115 5-20 | 5-21 | 32-03 
116 259 | 265 | | 43-73 
117 17 | | | 54-21 
4 118 | 0-58 | 0-58 69-57 
119 O83 | | 84-74 3 
3 
= 
a 


VOL. 
6 


956/57 


The quaternary system: benzene —carbon tetrachloride acetic acid water at 25°C 


The tie lines were located with the aid of two 
charts. The non-aqueous phase specific volume 
data were plotted on a triangular composition 
chart, vertices benzene, carbon tetrachloride, and 
acetic acid + water, 100cm sides, as lines of 
constant specific volume at 0-005 ml/g spacing 
over regions of interest. Specific volumes in the 
benzene — carbon tetrachloride — water ternary 
were calculated from the data of Scarcuanp et al. 
[19] on the densities of benzene — carbon tetrach- 
loride mixtures, solubilities in the system being 
estimated by interpolation between the binary 
solubilities. For the aqueous phase specific 


volume (1 cm to 0-002 ml/g) was plotted against 
weight per cent water (lcm to 2%). 

The phase compositions were determined as 
follows : 

(a) For all ternary tie lines the non-aqueous 
phase compositions were determined from the 
specific volume and isotherm data, and for quater- 
nary tie lines numbers 18, 19, 26, 27, 35 and 36 
from the specific volume, acetic acid concentration, 
and isotherm data. Knowing the composition 
of the initial mixture and that of the non-aqueous 
phase, the aqueous phase composition was 
obtained by trial-and-error mass balances. 


Table 2. Experimental tie line data for the ternary systems. Benzene (B)-— Acetic Acid (A) — Water 


(W) and Carbon Tetrachloride (C) 


Benzene Acetic Acid Water 


Acetic Acid (A) — Water (W) 


Initial mixture 


Non-aqueous phase Aqueous phase 


— Wt. aqueous phase 
WE. both phases Specific Specific 


Ww 
we % 
39-86 0-35378 
47°25 0-6400 
40-10 0-6393 

0-6381 
0-6652 


0-6601 


Volume Volume 


a25°C | 3 at 25°C 


1-14085 | 097108 


11414] . 0-97102 
113810 ‘96110 
1-13723 ‘95975 
110571 - ‘96354 
26-02 | 67-91 


110155 | 0-96611 
25-98 66-95 


Carbon Tetrachloride Acetic Water 


Initial mixture 


Wt. aqueous phase 
Wt. both phases 


O-4199 
O-4595 
0-4990 
0-5291 
0-4992 
0-5378 
0-5098 


0O-4885 


Non-aqueous phase Aqueous phase 


Specific Specific 
Volume Volume 


at 25 at 25°C 


0-63253 0-98226 
0-63479 0-96814 
0-63895 

0-64229 O-94576 

OG5617 O-92591 
0-607 67 0-91243 
0-67720 | O-90250 


0-69040 O-B8S54 


| 
Tie | --- 
line 
No. B A A 
2 35-39 17-36 
Rete 3 | 34-76 25-14 
4 34-72 26-56 
5 31-67 | 51-53 
6 32-13 52-21 15-67 | 
| 
> | | 
| 
| | 
No, | ( A A 
| 
7 6-42 6-78 36-79 - 
5 53-73 12-92 33-35 | 
10 Mi-47 26-87 26-66 
11 1-31 30°37 20-32 | 
12 | 39-53 1-92 19-56 
13 | “0-22 | 13-02 ~ 
144 67-37 
| SEIS | 40-88 7-07 
253 
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Table 3. Experimental tie line data for the quaternary system. Benzene (B)-— Carbon Tetrachloride 
(C) — Acetic Acid (A) — Water (W) 


Non-aqueous phase | 


Initial mixture Aqueous phase 


— Wt. aqueous phase 


Specific 


18 36-12 | 11-95 14-51 37-42 
11-61 22-33 31-03 
20 30-56 9-97 45-55 13-92 


| 
21 20-62 | 20-57 | 16-00 42-81 
4 22 20-06 | 20-04 24:52 | 35-38 
23 | 23-15 23-00 | 32-64 21-21 
M 24 17-68 17-62 | 48-7: 15-97 
= 25 | 17-97 17-97 | 51-25 12-81 
| | 
7 26 | 25-65 | 25-56 | 1365 | 35-14 
| | 
: 27 | 24-82 24-70 | 21-21 | 20-26 
3 
28 21-15 | 20-97 | 44-35 13-53 
4 29 20-80 | 20-62 | 47-81 10-77 
4 30 | 13-40 | 1484 81-69 
12-13 36-06 | 21-60 | 30-21 
32 12-53 8744 80-71 | 19-32 
33 10-11 30-11 | 4560 | 1418 
: OST 29-71 48-04 11-48 
; 35 13-78 41-11 12-84 | 82-27 
36 | 1-21 39-72 19-99 «| «27-08 
az | 1-77 | | 4163 | 11-60 
38 | 131 33-91 44-06 9-62 


0-5831 
0-5890 
0O-5124 
00-6476 


00-6562 


04894 


05516 


0-5666 
04583 
0-082 
0-4803 
0-6005 
0-6282 


O-4446 


0-4584 


0-5123 


00-5566 


1-01659 


1-01585 


1-00875 


0-88988 
0-89138 
0-89613 
0-89984 


| 0-90359 


0-89035 


0-89241 


0-90282 


0-90761 


0- 76401 
0-76693 
O-77467 
0O-789045 
0-79899 


| 076347 
| 0-76630 
0-79313 


0-80417 


0-97031 
2-36 
4:54 
0-95879 
483 
21-65 
0-95335 
21-44 


0-97007 
| O95773 
| 0-94635 

094122 | 
| 0-93968 


190 | 096072 
| | 
3-95 
18-21 

0-94147 


0-94022 


| 0-96521 
| 0-95618 
9-94154 
0-92704 


0-91826 


1-51 
0-96866 
1-47 
3-05 
0-95584 
3-05 
15-85 | 9.92580 
15-84 | 
20-91 
| 
20-96 


Tie Wee. both phases: Specific | 

line B Cc | A Ww Volume A Volume A 

No. | Wy | wey | we % | We % at 25°C wt. % | at 25°C wt % 

16 | 27-96 9-32 25-54 37-18 | 0-6141 | 1-01497 0-95890 | 

17 26-82 8-90 50-69 13-59 | 0-6279 1-00410 | 0-95629 


26-36 
26°35 
39-08 
39-10 
66-53 
66-86 


26-49 
26-32 
39°37 
39-40 
66-70 
67-06 
66-54 
66-59 


27-11 

27-10 
40-19 
40-13 
67-78 
68-08 
68-98 
67-60 


(b) 


ratio, by trial-and-error calculations. 


in Tables 4 and 5. 


The bi-nodial compositions obtained are shown 


For all other quaternary tie lines, the 
phase compositions were found from the non- 
aqueous specific volume, the composition of the 
initial mixture, and the observed phase weight 
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COMPARISON WITH 


PUBLISHED 


Benzene — acetic acid — water 
Hanp [9] has published careful and extensive 
data on this system. His saturation results show 
somewhat lower solubilities, but the difference is 
small, and of the same order as his random error. 
HaNv’s tie line correlation line is shown in Fig. 5, 


DATA 


q 
| 
0-5146 
| | | 
| VOL 
6 
| 1956 
| 
24-67 
| 
| | 
| 
| | | | 


VOL, 
6 


956/57 


The quaternary system: benzene -carbon tetrachloride - acetic acid —- water at 25°C 


Table 4. Equilibrium tie line compositions for the ternary systems. Benzene (B)-— Acetic Acid (A) - 
Water (W) and Carbon Tetrachloride (C)— Acetic Acid (A)— Water (W) 


Benzene — Acetic Acid — Water 


Non-aqueous phase Aqueous phase 


B A Ww 


No. | | Wey | | | 
| 
1 | 9722 | 266 | 0-12 057 | 25:56 | 73-87 
2 | 9720 | 268 O12 0-57 | 25-62 | 73-81 
: 94-84 499 | O17 | 0-93 36-49 | 62-58 
4 9420 | 561 | O19 | 106 | 38-31 | 60-65 
5 7477 | 2405 | 1-18 10-00 65-35 | 24-65 
6 72-67 | 2595 | 188 | 1124 | 65:73 | 23-03 


Carbon Tetrachloride — Acetic Acid — Water 


Non-aqueous phase Aqueous phase 


Tie | | 


Line 4 A Ww Cc A | Ww 
No. wey | wey | | | 
7 | 99-62 0-36 0-02 0-30 | 1512 | 8458 
8 | 9902 | 096 | 0-02 051 | 2697 | 72-52 
9 | 9793 | 204 003 40-91 58-02 
10 97-02 2-04 0-04 1-80 48-02 50-18 
11 95-44 | 449 | 007 | 3-24 56-22 | 40-54 
12 | 93-14 6-76 | 0-10 5-60 | 62-53 | 31-87 
13 89-90 992 | o18 | 9-50 | 66-39 | 2411 
14 87-21 12-54 0-25 | 12-77 | 67-27 19-96 
83-32 16-26 0-42 17-41 | 66-70 


Wt. fr acetic acid in aqueous phase 
Wt. fr water in aqueous phase 


+—- 


+0 0-0 
Wt. fr acetic acid in non-aqueous phase 
log Wt. fr benzene in non-aqueous phase 


log 


Fic. 5. © This work @ Bexrovrov 
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Table 5. Equilibrium tie line compositions for the quaternary system. Benzene (B)— Carbon 
Tetrachloride (C)— Acetic Acid (A)-— Water (W). 


Non-aqueous phase Aqueous phase 


Tie 
Line | B A Ww | B A | Ww 
No. | we.% | | Wey | | | % | 


16 | 23-87 487 OF15 | O84 | O18 | 38:53 | 60-45 
17 | 54-79 19-21 24-77 10-25 | 2-79 66-05 20-91 
18 73-25 24-26 239 | O10 O44 | O12 26-15 73-29 . 
19 | 71-31 23-71 O-15 0-84 0-20 38-84 60-12 
20 19-78 0-90 8-03 1-99 24-51 : 


21 49-00 49-05 1-89 0-07 O33 | O21 26-09 73-37 
22 47-90 48-20 3-79 O11 0-63 0-39 38-08 60-00 
23 45-95 8-57 O23 2-12 1-16 55-54 $1-18 
24 30-44 42-23 17-71 0-62 5-84 4-23 65-61 24-32 
25 36°45 39-52 23-00 1-03 8-29 6-68 66-05 18-08 


49-08 $8-05 0-07 O25 73-08 
47-06 47-99 3-04 O11 0-66 39-24 59-71 
39 47 0-65 $30 24-00 
36-35 38°27 24-23 1-15 8-01 7-12 65-85 18-12 


" 30 24-57 73-61 1-78 0-05 0-20 0-43 30-28 69-09 
q 31 24-22 72-70 3-01 0-07 0-43 0-60 39-59 59-38 
22-92 70-02 6-92 1-29 2-17 56-46 40-08 
4 33 20-10 64-64 14-85 0-41 3-46 7-14 66-06 23-34 
18-58 61-01 19-71 0-70 4-72 11-18 66-24 17-86 
35 24-67 73-79 1-49 0-05 O19 O35 27-40) 72-46 
4 24-08 72-80 3-05 0-07 0-38 0-69 39-08 58-95 
Be 20-19 3-63 0-45 66-29 2 
19-00 0-63 0-77 66-01 


together with the points representing the present 
q work. The agreement shown is possible, despite os 
the difference in isotherms, as the AND correla- | 
tion only involves the concentrations of the major ap 
components. Data of Bexrourov [2] are also 4 
shown ; his points lie on a line parallel to, but 3)20-0 7 
distinct from, the present results, and in view cle . 
of the coincidence of the latter and of Hanp’s sic . 
results, BEKTOURKOV’s are considered less reliable. rae 
che 
Carbon tetrachloride — acetic acid — water “= | 
> 
Smirh and Berman [23] determined this 
isotherm as part of a wider investigation, and 
their results are consistent with the present, - WL. fr. acetic acid in non-aqueous phase : 
except for the aqueous phase with acetic acid " Wt. fr carbon tet. in non-aqueous phase 7 
F concentrations below 50°, where their results Fic. 6. : 
appear subject to rather large experimental error. This work @ Bexrovurov A SALETORE et al 
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Saturation results of SaLETORE, MENE and 
WARHADPANDE [18] at 27°C scatter about the 
present data. They also reported three tie lines, 
shown on a Hanp plot together with the present 
results in Fig. 6. These lines were also determined 
by the synthetic method, and as the composition 
of one phase only of one tie line fell within their 
experimental isotherm, they must be subject 
to considerable error. 

Bektourov’s data [2] are also shown in Fig. 6, 
and again diverge from the present results. 

and Linke [21] quote saturation 
compositions ascribed to SicGi1a and Hanna [22]. 
These lie well within the homogeneous region 
as defined by the present results, and reference 
to the original publication showed that these 
compositions had no connection with the isotherm, 
and merely represented certain synthetic mixtures. 


CARBON 
WATER 


THe System BENZENE - 
TETRACHLORIDE 
The quaternary system studied includes three 
ternary systems with miscibility gaps : benzene - 
acetic acid — water, carbon tetrachloride — acetic 
acid — water, and benzene — carbon tetrachloride 
— water. While experimental results are presented 
for the first two of these, no attempt was made to 
obtain any for the last. This is a double non- 
consolute pair system with very low solubilities 
(the highest binary solubility being that of 
benzene in water, 0-18g per 100g), so that 
reliable results would be difficult to obtain, and 
it was not expected that they would be required 
in the present study. For many purposes the tie 
lines in this system may be regarded as passing 
through the water apex, and for cases in which this 
concept is inadequate, the distribution to be 
expected may be derived [17]. 


DiIscussION 

The results presented for the isotherm of the 
quaternary system show that the cloud point 
titrations were generally reproducible to within 
0-1% or better, and with the technique described 
it may be possible to improve the reproducibility 
even further by more careful work, e.g. repeat 
titrations. Such extra precautions were not 
taken in the present work as the nature of the 


benzene — carbon tetrachloride 
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-acetic acid— water at 25°C 
isotherm made correction of titration results 
possible by the methods described. 

Such precision is not, however, attainable over 
all the surface. Table 1 shows that for certain 
points larger corrections had to be made. These 
points lie in the region of the plait points, and 
during titration here unexpected phenomena were 
observed, which will be described briefly, as they 
As the titration 
proceeded a blue haze developed in the solution, 
which became marked, and rendered the solution 
almost opaque ; this was followed by the sudden 
appearance of a very heavy white cloud. As the 
development of the blue colour was very gradual, 
while the advent of the white cloud was most 
marked, and seemed reproducible, the latter 
event was taken to be the attainment of satura- 
tion. When, after weighing and slight heating 
to clear them, the solutions had been returned 
to the thermostat, they became definitely cloudy, 
but gentle twirling of the flask would be sufficient 
to clear them, suggesting that if the correct 
end point had been exceeded, it had been exceeded 
by very little. Plotting the results as above 
suggested that the over-titration was appreciable, 
and attempts were made to avoid this, by stopping 
the titration, and shaking, at frequent intervals 
after the first appearance of the blue haze, and 
awaiting the appearance of a definite cloud point. 
This occurred occasionally, but on the whole 
this technique was no more successful than the 
original. 


may lead to confusing results. 


As these phenomena were confined to 
the cloud point region, and as this region was of 
little interest for the purposes of the present work, 
no further investigation was made. 

Other workers, e.g. HaNnp [9], have drawn atten- 
tion to anomalous behaviour in this region, and as 
he pointed out, this emphasises the assumption 
implicit in the synthetic method that the cloud 
point curve is identical with the equilibrium 
saturation curve. This assumption is unnecessary 
in the analytical method, but that method can 
not often be applied, and generally lacks precision. 
This criticism would seem to be borne oui by 
the present results for the tie lines, where there 
are comparatively great discrepancies between 
the acetic acid concentrations, as found by 


titration, and those calculated. It seems likely 


) 
VOL. ; 
7 
956/5 
| 
| 


that the titrated, rather than the calculated 
figures were less accurate, as 


(i) The experimentally determined aqueous 
phase density checks with the calculated com- 
positions ; 

(ii) At the higher acid concentrations the 
repeat titrations are inconsistent, and several 
figures are above the limits of heterogeneity as 
found by the cloud point titrations ; 


(iii) A numerical check is available for those 
lines of group (a) above, for which a neat separa- 
tion was achieved, and for which the other data 
required are known to be reliable. For such 
lines the ratio weight of aqueous phase to weight 
of both phases was calculated from the phase 
compositions and the composition of the initial 
mixture. These figures, together with those 
obtained experimentally, are shown in Table 6. 
The difference between the two sets is small, 
and of the order to be expected from the methods 
employed. For eight of the eleven tie lines the 
calculated ratios are greater than the experi- 
mental values, but a statistical test [8] showed 
that there was no consistent difference, and the 
two methods of calculation should lead to 
consistent results. 

Extensive experiments were carried out to 
trace the source of the discrepancies in the acid 
concentrations. These showed that titration 
gave the correct amount of acid in the samples 


Table 6. Comparison of calculated and observed values of ratio 


R. G. H. Parnce and T. G. Hunter 


submitted, and that the overall composition 
of the equilibrium mixture remained constant 
during a run, so that the rise in acid concentration 
took place during sampling, by evaporation of 
the more volatile components, as the phase was 
sampled directly into a large vessel suitable for 
the titration. Avoidance of such an error would 
be difficult, and this would tend to restrict the 
precision of an analytical method. Although 
the acid titration errors were large compared 
to those of the methods adopted, they represent 
at most 8%, errors in the acid mass balances, 
and hence would have been considered accep- 
table in most liquid — liquid equilibrium work. 
In the method of tie line determination 
adopted, a sample was only required for the 
density measurements, and as this was pushed 
directly into a pycnometer, there was no similar 
opportunity for evaporation and consequent 
error. For the purposes of Table 6, acid titration 
results were accepted for acid concentrations 
below 5°, as here comparatively large samples 
were taken, and the errors involved were of the 
same order as those of the other measurements. 
Some estimate of the precision of this method 
of tie line determination may be obtained from 
Table 6. The average difference in the two sets 
of weight ratios is 0.00015, which corresponds to 
an error of 0-08°%, in the mass balance of a com- 
ponent in terms of the total weight. While this 
is certainly an underestimate of the average error 


Wt. aqueous phase 
Wt. both phases 


Wt. aqueous phase 


- calculated 
Line | Wt. both phases , 
2 | 0-6397 
+ 0-6385 
5 00-6654 
00-4597 


0-4991 
15 | 04881 
18 0-5100 
9 0-5147 
a | 0-4892 


35 | 0-4450 
0-4587 


Wt. aqueous phase 
Wt. both phases 


(observ. d) 


0-6400 

0-6381 
0-6652 
0-4595 
0-4990 
0-4883 
0-5094 
0-5146 
0-4894 
0-4446 
0-4584 


. 
- 
6 
1951 
| 
4 
258 


VOL. 
6 
1956/57 


The quaternary system : 


for all the tie lines (as those included in the 
table were necessarily not a representative sample) 
the average error should be of the order of 0-1°%, 
the figure originally aimed at. 

At the same time, the practical convenience of 
the technique described should be stressed, 
particularly for ternary systems, where, if the 
phase weight ratios are accepted as reliable, 
no samples are needed, and the only measurements 
required are the weights of the components added 
and of one phase. It is therefore suggested that 
this method represents advances over the con- 
ventional methods, both in convenience and in 
accuracy. 


RELATIONSHIP BETWEEN THE TERNARY 
AND QUATERNARY SYSTEMS 
BraNnckeR, Hunter and Nasu [4] have 
studied the quaternary system acetic acid — 
chloroform — acetone — water at 25°C, which is of 
the type illustrated in Fig. 1 (a). The quaternary 
tie lines were found to lie on the intersection of 
a plane passing through a tie line in a ternary 
system containing the pair of non-consolute 
components and the opposite apex of the tetra- 
hedron with a similar tie line plane based on the 
other ternary system containing the pair of non- 

consolute components. 

Pratr and Giover [16] determined tie lines 
in a similar system, actone — acetaldehyde 
vinyl acetate — water at 20°C, which seem to be 
in fair agreement with those expected from the 
Brancker, Hunter and Nasu _ relationship. 
BranckeEr [5] claimed to obtain better agreement 
by a different method of calculation, but this 
does not appear to be valid. 

No other quaternary systems of this type 
have as yet been investigated in sufficient detail. 
The application of the Brancker, HUNTER, 
and Nasu correlation to the present system gave 
systematic between actual and 
calculated results outside the expected limits 
of experimental error and it would appear that 
these tie line relationships do not apply to this 
system [17]. 

CHANG and Mouton studied four quaternary 
systems of the type illustrated in Fig. 1 (b). 
Of these the system water —ethyl isovalerate 


deviations 


benzene carbon 
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tetrachloride — acetic acid — water at 25°C 


ethy! alcohol — benzene was studied in some detail 
and its quaternary tie lines were found to lie in 
planes passing through tie lines in the ternary 
system containing the two non-consolute pairs. 
In Fig. 7 the line FG is a tie line in the two non- 
consolute pair ternary BCD, and the line EC 
is this line extended. The plane through EF GC 
and the opposite apex A, namely plane ECA, 
is then a quaternary tie line plane, i.e. a quater- 
nary tie line through any point in the plane lines 
within that plane. Cuanc and also 
obtained distribution data for three similar 
systems, which supported the existence of such 
quaternary tie line planes, within experimental 
error. 


Prince [17] recalculated Brown's [6] quater- 
nary tie line data for the system aniline - 
— cetane 


benzene 
cyclohexane at 25°C and concluded 
that the tie line relationship suggested by Cuanc 
and Movu.tton might apply to this system. 
Moreno and GonzaLes [12] investigated the 
system methanol — carbon disulphide — oleic acid — 
olive oil at 25°C and treated this as if it were a 
quaternary system of type (b). They claimed that 
the CuanGc and Mou tron relationship applied to 
this case. 

If the Cuanc and Moutton relationship applies 
to the present system then the quaternary tie 
lines must be in planes through the water and 
acetic acid apices, namely planes of constant 


| | 
¥ 
Fig. 7. 
de. 
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Table 7. 


R. G. H. Paice and T. G. Hunter 


Ratio carbon tetrachloride to benzene for quaternary tie lines 


Tie Overall 
Mixture 


18 00-3308 
19 ©3315 
16 0-33835 
20 0-3262 
7 0-3318 
21 0-9976 
26 0-9965 
22 
27 00-9952 
23 0-9935 
24 0-9966 
28 0-0015 
25 1-0000 
29 
35 2-983 
30 2-000 
31 2-973 
36 3-007 
32 2-988 
33 2-978 
37 2-074 
34 3-010 
2-946 


Non-aqueous 


Aqueous phase 


Phase 


0-3312 0-278 
0-3325 0-238 
0-3357 0-214 
0-3396 0-248 
00-3506 0-272 
1-0010 0-636 
0-9973 0-833 
1-0063 0-619 
1-0006 0-591 
10155 0-547 
10707 0-724 
1-0435 0-713 
1-0842 0-806 
1-0528 0-799 
2-991 1-84 
2-996 2-15 
3-002 1-39 
3-023 1-82 
3-055 1-68 
3-216 2-06 
3-152 2-07 
3-284 2-37 
3-138 207 


carbon tetrachloride ‘benzene ratio and _ this 
ratio must be the same for both conjugate phases 
and for the overall mixture from which the 
conjugate phases were derived. If in Fig. 7 
the apices A, B, C and D represent acetic acid, 
benzene, water, and carbon tetrachloride, res- 
pectively, then the Cuanc and MouLtTon quater- 
nary tie line plane will be the plane ACGFE, 
and in this plane the ratio of B to D, namely 
benzene to carbon tetrachloride, will remain 
constant. 

In Table 7 this ratio, as carbon tetrachloride to 
benzene, is given for the overall mixture and for 
the conjugate phases of the quaternary tie 
lines, arranged in order of increasing acetic acid 
content. This ratio deviates in both phases in a 
systematic manner from the overall ratio, while 
there is no such deviation in the results of CHANG 
and Mouton. As benzene and carbon tetrach- 
loride have very different densities (0-87 and 
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1.59 respectively) and as density was employed to 
determine the ratio of these components in the 
non-aqueous phase, the determined ratio data 
for this phase are believed accurate. In the 
aqueous phase, however, the concentrations of 
both these components are low, and the ratio 
of these in this phase will be subject to fairly 
large errors. 

As far as the present system is concerned it 
would appear that the Cuanc and Mouton 
tie line correlation is only an approximation. 
Even if this, or a similar correlation were valid 
generally for quaternary systems of type 1 (b), 
it would not be sufficient to predict quaternary 
equilibria from data for the ternary systems 
only (as is possible with the Brancker, HUNTER, 
and Nasu correlation for type 1 (a) systems). 
Correlations which enable such predictions to be 
made have been discussed in more detail by 
Prince [17}. 
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On shape factors for irregular particles—lI 
The steady state problem. Diffusion and reaction 
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Abstract—The problem of heat or matter transfer for a bed of spherical particles has been 
fully solved by Amunpson. In this paper we consider the modifications which must be made in 


steady state solution when the bed may consist of irregularly shaped particles. Considering 
the case of diffusion limitation of a first order reaction it is shown that the results for all shapes 
will lie close together if the characteristic dimension of the particle is taken to be Up/ Se: the ratio 


of its volume to its external surface area. 


Résumé—Le probléme de transfert de chaleur ou de matiére pour un lit de particules sphériques 


a été résolu par AMUNDSON. Dans cet article, 


lauteur considére les modifications susce ptibles 


détre apportées a la détermination ce l'état de régime quand le lit est constitué par des particules 


de formes irréguliéres. En considérant le cas de limitation par diffusion d'une réaction du 1° ordre 


auteur montre que les résultats sont voisins, si la dimension caractéristique de la particule est 


INTRODUCTION 


Tue notion of a shape factor for particles of 
irregular shape has been introduced empirically 
in several contexts interesting to the chemical 
engineer. Leva [6], for example, in studying 
pressure drop in packed beds introduces the shape 
factor ¢ — (36 7)’ v,? s, where v, is the volume 
of the particle and s, its surface area. This factor 
is unity for the sphere and between 0 and 1 for 
ali other shapes, and using it he is able to improve 
the correlation of pressure drop versus ReyNoLD’s 
number. The introduction of such a factor should 
be somewhat tentative for mathematical con- 
siderations give no hope of a single parameter 
completely reconciling the solutions of a particular 
problem for two distinct shapes of particle. 
However the practical value of so simple an 
adjustment of results is very great even though 
a small margin of error is introduced. It is our 
purpose here to examine two problems in packed 
beds and show the extent to which such a practice 
can be justified. 

The first problem concerns the diffusion limita- 
tion of a first order reaction taking place on the 
surface of a porous catalyst. This problem was 


égale a f S, (rapport de son volume a sa surface externe). 
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first solved by Ture te [8] who gave two curves 
one for flat plates and another for spherical 
particles. His solution and the physical back- 
ground has received an extended treatment by 
Wee er [9] and is described and used in HouGEN 
and Watson's well-known treatise [4]. In sections 
2, 3 some justification is given for using as the 
characteristic linear dimension of an arbitrary 
when a single 
curve for the effectiveness factor may be used 


shaped particle the ratio 


with good accuracy. In section 4 the equivalent 
particle diameter for a mixture of particles is 
discussed. 

The second problem, to be considered in Part 
II, is that of the fixed bed heat exchanger or ion 


exchange column. It shows well the pitfalls of 


endeavouring to use a single shape factor, for 


though at first sight the same specification of 


characteristic length appears to give a measure 
of agreement between various shapes, this 
agreement is shown to be spurious. However 
a shape factor can be found in this case which 
compares the * softening " ofa sharp temperature 
or concentration wave passing through the bed 


for various shapes of particle. Since a very 
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general problem has been investigated by 
AmUNDSON [1] for spheres we follow his notation 
Section 
7 gives a method of calculating the shape factor 


for this problem and the final paragraph discusses 


and methods some distance in section 6. 


the relation between the two problems treated 
here and that solved completely by AmuNDsON 
for spherical particles. 

DIFFUSION AND REACTION 
IN CATALYST PLATES 


2. SIMULTANEOUS 


To introduce this problem we will first give 
Tuie.e’s solution for the flat plate in terms of 
the notation to be used here. Consider a porous 
catalyst particle in the form of a thin flat plate 
with sealed edges so that the flow of reactants is 
Take a co- 


ordinate x in this direction of flow and let z = 0 


always perpendicular to the faces. 


denote the centre plane of the plate and 2 L @ 
the faces. 


order reaction with no volume change whose 


We consider for simplicity a_ first 


rate is proportional to the concentration c¢ of 
one of the reactants. Let 


c ¢ (av) concentration of the reactant in 
the plane 
e € (a) its concentration at the surface 


D — diffusion coeflicient of the reactant 
within the particle 

o = active catalyst surface unit 
catalyst volume 

k reaction rate constant unit active 


catalyst area. 


Taking a matter balance over an element between 


the planes 4 } dx we derive in the usual way 


d*c 
) 
da? 


koc (1) 


By symmetry there is no flow of reactant across 
the centre plane 


a 0 (2) 


and at the surface 


These equations determine the distribution 
of concentration within the particle and we have 
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A 
where A =ko/D (4) 


s 
cosh Aa 


and the flow into the particle per unit area is 


|p | Dc, A tanh Aa. (5) 
dz lene 


Consider now a bed of such particles through 
which the reactants are flowing, the particles 
being sufficiently small for the concentration c, 
to be constant around any single particle. 


Let z = co ordinate of bed depth 
A = area of bed 
x’ = voidage of bed 


vy = volume rate of flow of. reactants. 
Then a matter balance over an elementary 
section of the bed between the plane z + } dz, 
gives 
[v, + 


a 


[vy 4 dz 


D 


In deriving this we have made use of the facts : 


(1 — «') Adz = volume occupied by catalyst 
particles 
a’ area through which reactant 


flows in ‘volume of catalyst. 


Using equation (5) and taking the limit dz + 0 
we have for c, the equation 
uy de, l dD 

(47 A tanh Aa) ¢,. 6 

A(l x') dz a (6) 

If there were no diffusion limitation and all the 

catalyst surface area were available to the reac- 
tant a similar matter balance would give 

Uy de, 


All — a) dz 


Comparing equations (6) and (7) we see that full 
allowance can be made for the diminished rate 


of reaction if in (7) & is replaced by nk, where 


tanh Aa 


8 
ia (8) 


7 is called the effectiveness factor. 


> 
| 
| 
4 
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Cast in this form the problem is immediately 
amenable to generalisation. We may set up the 
equations for a particle of arbitrary shape, to 
replace equations (1), (2 and (3) for the flat 
plate and then calculate the average rate of flow 
across the particle surface to replace equation (5). 


REACTION IN CATALYST 


ARBITRARY SHAPI 


3. AND 
PARTICLES OF 
Consider a catalyst particle of arbitrary shape, 
but as before sufficiently small for the concentra- 
tion of the reactant, c,, to be sensibly constant 
over its surface. Let &, », 
ordinates in the space of the particle which 


{ be Cartesian co- 


we take to occupy a region R and be bounded 
by a surface S. These space variables correspond 


to the variable x above. 


Fic. 1. The general particle 


Let v, volume of particle, 
di element of volume, 
s, = external surface area of particle, 
ds element of surface area, 
) = differentiation along the outward 


normal n, to S. 
v? the Laplacian operator, 


A matter balance within the volume of the 


particle now gives for the concentration 


¢) the equation 
Vic = A*c in R, (9) 


where again A* =koe D, and on the surface 


c c, on (10) 


Equations (9) and (10) replace (1)-(3) and their 


solution gives the concentration distribution 
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within the particle. The total flow into the particle 


over the surface is 


{ 


R 


D | ds ko | | edv (11) 

s R 
by Green's formula and the equation (9). 

Let x (€, ». ©) denote the solution of (9) which 
has boundary values xX 1 on SS. Then 
e(é,», 0) =e, xX (é 9,0) and a matter balance 
similar to equation (6) gives 

uy de 
A (1 x) dz 


| X to| koe, (12) 


from which it is evident that 
irre 

| | x dv (18) 
“pede 


MxinR 


lon 8S. (14) 


where 


and x 


The effectiveness factor » depends on the shape 
of the region R and on the parameter A a, where 
a is some characteristic dimension of the particle. 
The problem of a shape factor is the problem of 
choosing the value of a appropriate to the given 
shape of the particle so that the functions » (Aa) 
shall be as close together as possible for all shapes. 

For certain common shapes equations (13) and 
(14) can be readily solved, taking for a an obvious 


dimension. For example 


(i) Flat plates: a = half thickness : 


» = tanh Aa Aa 


(ii) Cylinder with ends sealed : a = radius : 
” 21, (Aa) \a I, (Aa) 
(iii) Sphere: a radius : 
3 (Aa coth Aa 1) A? 
(iv) Cylinder with porous ends: a radius, 
length = 2a p: 
l 
1-- ) >" 
7 (2n + 
Va? + j2+(2n4+ 1) p* x 


In the above J, (x) is the modified Bessel function 


: 
a 

6 
956 
| Vv? x 195¢ 
\ 
4 
4 
ne, 
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of order n, and j,, is the m™ zero of the Bessel 
function J, (2). 

When Aa + 0, » +1 as it should, for Aa = 0 
implies no limitation by diffusion. When 
Aa >a, » >m/Aa where in the first three 
cases we have m = 1, 2, 3, respectively. Now 
in all three cases a/m is the ratio v,/s, so that for 
large Aa, » ~ (A v,/s,)'. Thus in three important 
cases we see that by setting a =v,/s, the 
functions » become identical for large and small 
Aa. Before suggesting that it is satisfactory always 
to set a = v,/s, we should do two things. Firstly 
estimate the error likely to be introduced for 
intermediate values of Aa and secondly show that 
for an arbitrary shape we always have 7» ~ 
(A v,/s,)". 

For the first let us set 


A = Xv, /s, = (v, ko/s2 
then for 


(i) flat plate: » = tanhA/A (15) 
(ii) eylinder: =I, (2A)/A1, (2A) (16) 


(iii) sphere : = (8Acoth38A — 1)/8 A? (17) 


Table 1. Effectiveness factor » for various shapes 
A 
Flat pla’e Cylinder Sphere 

0-1 997 995 904 
0-2 ‘O87 ‘981 ‘O77 
0-5 924 “892 ‘876 

762 698 O72 
2 482 432 | -416 
5 197 | 
10 100 100 97 


From this table of values it is clear that the 
functions 7 (A) follow one another quite closely, 
the greatest difference between spheres and flat 
plates being approximately 0-09. To obtain 
estimates of the divergence in the general case 
would be an excessively difficult task, equivalent 
to solving the problem completely, but this 
table gives us some confidence in the use of the 
parameter A. 

As an illustration of the variation of » with 
shape at constant A we may mention the spherical 
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1 — 
os — at plate] — 
Cylinder + -+4 
ot 4 Sphere 
| 
0-1 
— 
| 
0-01 2°0 10°0 50-0 
A 
Fic. 2. Effectiveness factor for various shapes 
shell. In this case we assume that the active 


catalyst is a spherical shell between the radii a 
andpa. Thens, = 4 andv, = 4 — p*)/8 
so that A = Aa(1 — p*)/8. Setting (1 — p*) A’ 
= 8(1 — p)A we have 


_B8(1—p) 1 
4? 
(1 — p)* (A’ — tanh 4’) + pA® tanh A’ 


(1— p) tanh A’ + pA’ 


When p =0, A’ = 3A and 7 reduces to the 
expression for a sphere, (17), whilst when p = 1 
we have A’ = A and the expression for a flat 
plate, (15). Taking A = 1 at which point we have 
noticed the greatest difference, we get the follow- 
ing values : 


p 0 25 ‘5 75 1-0 
672 ‘675 691 723 ‘762 


As regards the asymptotic behaviour of 7 it is 
at once evident that A = A v,/s, is likely to be a 
For large values of A 
imply a strong diffusion limitation, i.e. the reac- 


uniformising parameter. 


tion rate is so much faster than the diffusion rate 

that the concentration falls off extremely rapidly. 

In these circumstances when only a thin surface 

layer is really operative it is effectively a flat 
irre 


plate, and » 


J. 


| x dv, receiving its main 
k 
contribution from the immediate neighbourhood 


of the surface, becomes proportional to s,/v,. 


| 
he fi 
ate 
‘ 


More rigorously, we consider § to be composed 
of a finite number of smooth surfaces as, for 
example, is the finite cylinder with its curved 
surface and flat ends. In the neighbourhood of 
any smooth part of S we may introduce coor- 
dinates such that S is the surface of constant 
(say = a) then = dé and a solution 
valid in the immediate neighbourhood of S for 
large A is 
x = e~ Me 4 (é, 

where ¥ and ¥, are negligibly small on S. Since 
drops to zero very rapidly it is sufficient to have a 
solution valid in the immediate neighbourhood of 
S and 


” x dv = V? x dv 
Up Av, 
R R 
p (18) 
1 ax 1 
ds | ds 
Me, J J dé Av, 


These considerations give a fair measure of 
confidence in using the mean curve of figure 2 
for all shapes. In this figure the curves of table 1 
are given to give some measure of the possible 
error. 

A remark may be made as to the use of the 
shape factor ¢ (36 s, introduced by 
Leva. If the general particle of volume v, is 
replaced by a sphere of radius a of the same 
volume, i.e. a = (3v,/4 7)’ and the parameter 
A* Aa 3 for the sphere is used, no allowance 
has been made for the difference in shape between 
the particle and the sphere. If however we take 


A = $A* = A(80,/4 (36 = Av,/s, 
we arrive at our parameter. On a log-log plot 
of » versus A such as in Fig. 2 the points 
corresponding to A for the particle and A*, the 
value of A for the sphere of equal volume, are 
separated by a constant difference log ¢. It is 
clear from this that of all catalyst particles of 
given volume diffusion limitation is most severe 
in a spherical one. 


4. Mixtures or ParTIcCLe Size 


If the bed is composed of a mixture of particles 
of various sizes we may enquire whether it is 
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possible to use a mean particle size to describe 
the overall performance of the bed. Leva in his 
studies of pressure drops [6] suggests a mean 
' = Dw,a;' where w, is 


i size and a, the 


radius a, given by a 
the weight fraction of the 
radius of the sphere of equal volume. We shall 
see that this is satisfactory in so far as any mean 
can be satisfactory. The bed is assumed to be 
homogeneous in the sense that a sample from 
any part of it would show approximately the same 
proportion /; of the packed volume to be occupied 
by particles of the i™ size. We have Df; =1 
and if the density of all particles is the same 
=f;. 

To account correctly for the mixture of sizes 
we must take each size and calculate the corres- 
ponding value of A say A,, then 


= (Aj (19) 


If this value of » corresponds to a particular 
value of A, say A, we have from the relation 
v,/8, A/\ a value for the equivalent mean size 
of the particles. Equation (19) is however 
quite general and can be used for a mixture both 
of shapes and sizes, nor should it be excessively 
difficult to use in practice. 

There are however two cases in which a mean 
size can be calculated before finding the effec- 
tiveness. The first is the trivial case when 
diffusion is not limiting for any size of particle, 
A < 0-1 for even the largest particle. In this 
case any mean size will serve for » will be 1 in 
all cases. The second case is when diffusion is 
limiting even in the smallest particles, i.e. 
In this case »(A,;) = A;' 
and as we have seen A ¢é A* ¢Aa/3 where 
a is the radius of the sphere of equal volume. 
From (19) 


A;' 


A> 5 for all sizes. 


(3/A¢) Df, a;' = 8/Ada 


where a is the correctly defined mean value. 


Hence 
a =(Zf,a))" (20) 


is the appropriate mean. 

If the range of sizes involves that diffusion is 
limiting in some particles but not in others then 
any mean size calculated before the effectiveness 
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will be inaccurate. For example if the mixture is 
in equal volume proportions of five sizes, f, = 0-2, 
giving A; = }, 1, 2, 5, 10, the value of A calculated 
from the mean (20) is 1-316. This gives 7» = 0-658 
which is widely different from the value 0-397 
calculated from (19). 


respectively, a practice which equation (21) 
justifies in this case. 


CONCLUSION 
It has been shown that, if v,/s, is taken as the 
characteristic dimension of an irregularly shaped 
particle, the solutions of the steady state problem 
of simultaneous diffusion and reaction are largely 
independent of the particle shape. For the tran- 
sient problem it will be shown in Part II that this 


5. ExTerRNAL Dirrusion LIMITATIONS 


A remark may here be added, in view of its 


’ We assume that because of the presence of a C = concentration in fluid (or cy) 
V thin film (of thickness h) of stagnant fluid around c, = surface concentration 
956/57 the particle the surface concentration c, is not D = diffusion coefficient of reactant in catalyst 


connection with the problem which follows, 
as to the additive effect of a diffusion limitation 
through the relatively stagnant film surrounding 
a pellet in the fluid stream. This is often termed 
a mass transfer limitation and has been con- 
sidered by Houcen and Warson [4] and others. 


quite equal to the concentration c, in the body of 


the flowing reactant. If h is small we may assume Dy = diffusion coefficient of reactant in fluid 
that the flow of reactant is linear and if z is a B = overall cilectivences factor 
coommnate through the thickness of the film fg = volume fraction of ith sine 
D, 0, where D, is the diffusion coefficient h = film thickness 

Io, 1, = modified Bessel functions 


in the film. Then ¢c =¢, 27 
and ¢(x) = (#/h)(c,—¢;) + ¢%. Thus for the 
inward flow through the film we have 


= weight fracti 
A matter balance in the bed then gives er eae 
z = length coordinate 
Uy dey = length coordinate 
A(l— a’) dz a’ = voidage of bed 
where E, the overall effectiveness is given by = Caster 
Y= film effectiveness factor 
4 ake (21) A = (ko/D)s 
E 8,D, 7 A = 0, A/s, 

In such a situation as this the reciprocal of the p = fraction of particle radius 
effectiveness is often called a resistance and the active catalyst surface/unit volume. 
total resistance E™ is said to be the sum of the = shape factor 
resistances 7", n,' due to pore and film diffusion x = solution of (9) 
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characterisation is inadequate, and a new shape 
factor will be calculated. 


= area of bed 


= Bessel functions 


= reaction rate constant 


= particle volume 


NOMENCLATURE 


radius of sphere or cylinder 


particle 


root of Jy (x) = 0. 


ratio radius : length of cylinder 


region occupied by particle 


outer and inner surfaces of R 


external surface area of particle 


volume rate of fluid 


A 
| 
Jo. J, 
in 
— D, = D, = nkee, 
v, hke - 
Thus c,=c,{1 + and 
8,D, 
nkoctc,=keo + 2 _ Cy. 
| 
| 
4 
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Studies in filtration 
The initial stages of the cake filtration 


HEERTJES 


Laboratory of Chemical Engineering of the Technical University Delft 


( First received 11 January 1957; in final form 25 February 1957) 


Abstract—The specitic resistance of thin cakes has been investigated, This resistance (r) besides 


depending on the system is influenced by the initial filtration velocity and by the concentration 


(c) of the suspension, showing a decrease by increasing the velocity and the concentration. 


An equation for the relation between r and ¢ at constant initial filtration velocity has been 


viven. 


Résumé—L ‘auteur a étudié la résistance spécifique de gateaux minces. Cette résistance (r) par 


ailleurs dépendant du systéme, est influencée par la vitesse de filtration initiale et par la con 


centration (c) de la suspension. Elle décroit quand la vitesse et la concentration augmentent. 


INTRODUCTION 


Ir has been described in an earlier publication 
[1] that for the system puritine~—nylon cloth the 
type of filtration is governed by the concentra- 
tion of the suspension. Increasing the concen- 
tration (c) from a value of about Il mg_ per 
litre showed first a blocking filtration [2] and 
than a cake filtration with a decreasing specific 
resistance of the cake (r) which resistance reached 
a constant low value. For the systems puritine 
twillfill cloth and glass-spheres-nylon cloth no 
blocking filtration was found, but the influence 
of the concentration of the suspension on the cake 
resistance appeared to be the same as given 
above. 

The main purpose of the work to be described 
here was to evaluate an equation for the relation 
between r and ¢ and also to check for some other 


systems the validity of the observations made. 


EXPERIMENTAL 


(With R. de Borer, D. H. Jas and H. van Veen) 


The experimental techniques and procedures 
used were the same as described in the two 
other papers [1, 2] to which the reader is referred. 


In a single case a small wooden plate and frame 


Lauteur donne une équation pour la relation entre r et ca vitesse de filtration initiale constante 
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filterpress with three plates has also been used. 
A piston-pump transported the suspension in 
that case. 

The suspensions were in water, the temperature 
was in all cases 20°C. For every experiment the 
normal precautions to prevent the filter effect 
were taken. 

The substances used were: two types of 
puritine: Pur [ [1] and Pur IIL (10-300 ,»); 
Pur I A (75-175 p); 
Pur I A after dry elutrimtion for 1 hour: Pur 
1B; Pur I B after dry elutriation for another 
hour: Pur 1 C; the elutriation product obtained 

Pur I D (< 15 two sieve fractions 
of Pur II: Pur II A (75-175) and Pur II B 
(175-300 Norit (also an active carbon, 


a sieve fraction of Pur I: 


this way : 


containing ZnCl,) and three sieve fractions : 
Norit A (> 200, 27-5 per cent), Norit B 
(100-200 36-0 per cent) and Norit C (2-10 
36-5 per cent); glass pearls as used before [1] 
(4-40 x): Glass p.; sieved aluminium oxide 
(75-175 w) : Alum. 


*In the puritine sieve fractions a rather large amount 
of particles much smaller than the smallest diameter of 


the sieves used occurred. These could be removed by 


elutriation with air 


| 
ay 
VOL. 
956/5 
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Two different types of filtercloth have been 
used : nylon cloth (I) with 4000 openings per 
em® and twillfill cloth (II) for which no reliable 
number of openings can be given. In_ the 
turbulent regime of flow the resistances (see 
below) were 0°56 and 2-0 respectively. 

In the following the filters and cakes will be 
characterized by their resistance to flow (R, for 
the filter, r for the cake (gram per em*) and 
R for cake filter). As the conditions chosen 
for the filtration pressure and the resistance of the 
filter and the amount of cake were such that 
always turbulent flow existedt, the resistances 
R and R, have been defined and calculated, 
by means of the equation ; 


1dV (AP)? 


1 
Adé R (or Ry) (1) 
Moreover : 
A aR 
(2) 


Therefore r has been calculated from the slope 
of the straight R-V line. The viscosity terms 
have been included in the resistance terms for 
convenience sake. This could be done because 
the experiments have been carried out at the 
same temperature of 20°C. Correction for the 
viscosity (») can be made by dividing r and 
R, by 4°"'. The pressure difference will be given 
in atmospheres (atm), length in cm and time 


in sec. Therefore R (or Ry) is expressed in : 


sec cm sec 
_ and r in 


em 


RESULTS AND Discussion 


The results obtained and the discussion of 
these results will be presented in two paragraphs. 


The influence of the initial filtration velocity on r 

From earlier work indications were obtained 
that the cake resistance was not only influenced 
by the concentration of the suspension but also 

tit may be remarked that this also holds for the 
experiments of the first paper [1]. These have been 
calculated assuming laminar flow. The values given 
there should therefore be multiplied by (0-17°* = 2-8 
The conclusions drawn there remain the same. 


M. Heertyes 


by the condition of the filtercloth. If this cloth 
was not cleaned thoroughly, no reproducible 
results could be obtained. Because the building 
up of the cake and therefore its resistance, 
for not too thick cakes, formed in not too long 
a time (with thick cakes other phenomena may 
occur such as a changed packing by continued 
flow or by diminuation of the flow velocity), 
is governed by the packing of the initial layers 
of the cake, a series of experiments have been 
carried out in which the resistance of rather 
thin layers of cake formed under different 
filtration conditions has been measured. Either 
the filtration pressure was changed and the 
cloth cleaned thoroughly before every experi- 
ment, or a series was carried out in which the 
cloth was only cleaned superficially by rinsing 
it in streaming tap-water after each experiment, 
thus increasing the cloth resistance because 
ever and anon some particles will remain in the 
cloth. In a few cases changes in AP and 
in R, were combined. Some of the results are 
given in the Figs. 1, 2 and 3. In all cases r is 
presented as a function of the initial filtration 
velocity ; 


vy, — 
0 
800) 
700) 
| 200 
© S00) 
400}- 
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Fic. 1. The cake resistance r of Pur I and Pur II as a 
function of the initial filtration velocity Ug for different 
concentrations. 

Substance Filter Concentration. 
Purl 110% 


@ Purl I 
¢ Purl I 510° 
+ Purl I 
> Purl I 10-10° 


> Pur WU 
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in which R, is the resistance of the clean or 
incrustated filtercloth. It has been assumed that 
the relation between r and U, at a constant 
concentration can be expressed by : 


r=a-+ BU,. 


This relation has no physical basis at all and is 
used only to facilitate comparison. Straight 
lines based on this assumption, calculated by 


Fic. 2. The cake resistance r of Norit and of a mixture of 
Norit A and Norit C, as a function of the initial filtration 
velocity Up. 

Norit 


x = 96 per cent Norit A + 4 per cent Norit C 


The initial stages of the cake filtration 


the method of the least squares have been 
drawn in the figures and evidently only give an 
approximative correlation and only between the 
limits of U, as measured. In Fig. 1 the results 
for Pur I and Pur II are presented ; here only 
the filtration pressure has been changed. In 
Fig. 2 the results for Norit and a mixture of 
96 per cent Norit A and 4 per cent Norit C 
(weight basis) are presented; here only R, 
changed. In Fig. 8 the results for Pur I are 
given; here both the filtration pressure and R, 
have been altered. In Table 1 the calculated 
values of « and 8 have been assembled. 
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U 
Fic. 3. The cake resistance r of Puritine I at two different 
concentrations as a function of the initial filtration velocity 


Ue 


Pur I I 10-10° 


Applied 
change 
AP 140 4 
8-105 (range : 250 115 
| 510° 0-07 — 0°25) 190 + 100 
s-10°% 830 430 
1-10°5 65000 69000 
Pur II I 1-10°5 785 400 
Pur Il Il 1-105 1140 540 
Norit | I! | 10-10° 190 265 
96 per vent il 10-10°5 Ky 120 200 
Norit A (range : 
+ 4 per vent O-5 -> 1-1) 
Norit 
Pur I il 10-10°5 AP (range 415 380 
20-10°% 0-1 0-45) 190 110 


and Ry (range : 
0-56 —» 2-02) 


, 
150 
ry 
Yo 
| ° 
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Some experiments have been carried out in a 
filterpress with Pur I. It was not possible in 
this press to carry out a filtration at constant 
pressure. Therefore the press was operated 
under a more or less constant rate and the 
corresponding filtration pressure read. The varia- 
tion in the pressure at each stroke of the pump 
was rather large, because no buffering zone was 
used. Again r has been plotted versus a U5, 
which in this case has been taken as the velocity 
for the clean cloth corresponding to the filtration 
pressure measured. The constants calculated 


were : 
C= 510° x = 195 8 = + 81 
= 10-105 = 185 L 16 
- 20-10% 180 L 


Uj range = 0-7 + 2-0 
=0-5 +30 


It is apparent that the accuracy of several 
experiments is not very great, especially of those 
at low concentrations and with the artificial 
mixture. This is not only caused by inaccuracies 
in the measurements, which are especially large 
at high values of r, but also by the experimental 
procedure, in which the apparatus has to be 
filled with water before the actual experiments 
start. This gives rise to a small volume region 
where the concentration changes, which may be 
accompanied by a rather big change in r in the 
beginning of the filtrations. The conclusions 
which can nevertheless be drawn from the data 
presented are that r depends on the initial 
filtration velocity, therefore on the velocity 
with which the particles arrive at the filter cloth. 
In general a cake formed at a low velocity com- 
pared with a cake formed at a high velocity, 
will show a denser packing and will therefore have 
a higher resistance, as found. An exception forms 
the experiment with Pur I over the nylon filter 
at a concentration of 5-10°, where the reverse 
is true. It is remarkable that using this concen- 
tration the experiments with Pur I always 
turned out to behave irregularly. An explanation 


has not been found. 


P. M. Heertses 


\ second conclusion is, that the influence of 


the filtration velocity on the cake resistance 1s 
more pronounced if the cake resistance is large, 
viz. if the packing of the cake is dense. Therefore 
for cakes formed from suspensions of a relatively 
high concentration the influence of the filtration 
velocity on the resistance of the cake is small 
and even can disappear. 

The above is partly confirmed by the experi- 
ments in the filterpress with Pur I. 

In analysing these results it has to be kept in 
mind that the first layers of the three cakes are 
formed at the lowest U, given and that the 
influence of U, as presented is not directly 
comparable with that of the experiments described 
before. Moreover whether the cakes are really 
compressible or whether the pulsating action of 
the pump caused a consolidation of the cake, 
“annot be decided here. However this parallelity 
with the former experiments exists : the influence 
of the increasing pressure difference — small in 
itself —- diminishes with an increase of the con- 
centration of the slurry from which the cakes 
are formed that is, diminishes with a decrease in 
the density of the packing. 


The form of the r-c line 


From paragraph (a) it follows that a reliable 
r-c relation can only be found under such con- 
ditions that always the same initial filtration 
velocity is guaranteed, which experimentally is 
most easily realized by using the same filtration 
pressure and a clean cloth. For nylon it was easy 
to clean the cloth, twillfill cloth was not so easy 
to clean causing greater uncertainty and greater 
deviations in the results. Precautions in this 
respect have always been taken as good as 
possible. 

If by c, is denoted the concentration at which 
the blocking filtration changes to the cake 
filtration, it has been found earlier for puritine 
that the resistance of the cake at this concentra- 
tion (r,) was about equal to the differential 
specific resistance of the blocking filtrations for 
g =0. Therefore, including the results formed 
for blocking filtrations with puritine [2] : 


R, 


~ Wa) aN 0 
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Table 2. Filtration of Pur II over twillfill filtercloth at AP = 0-21 atm 


| 


No. Concentration |r (atm.°55 sec cm/g) 


(gram /litre) 


0-019 400 
0-039 214 
0-015 430 
0-024 230 
0-034 186 
0-042 104 


0-033 206 
8 0-043 184 
9 0-047 184 

10 0-050 180 


0-032 
0-020 
0-020 


208 
276 
280 


0-014 450 
15 0-010 860 
16 0-010 780 
17 O-O15 1260 


0-024 


OW) 


Ry (atm.°'55 sec /cm) Ry, Ry 


| 3-2 | 1-75 
24 18 
| 2-0 2-0 


2-3 | 2-1 
2-4 1-75 
+> (1) > 1-2 7-0 


From r, onwards the resistance r decreases 
rapidly with an increase of the concentration to a 
rather low value, which upon further increase of 
the concentration remains constant, that is 
deviates not more than the experimental error 
from this value. The density of the packing of a 
cake, and therefore its resistance, at a certain 
filtration velocity, will depend on the amount 
of particles arriving per pore in the same time 
interval. Therefore r will depend on ~ For 
the type of relation between these two magnitudes 
a probability distribution seems likely, which will 
be given in the form : 

r—r.=e 4, in which r_ is the value of 
the resistance at the concentration ¢c = co. As 
r =r, fore = ¢, the relation can be changed to : 


m an (3) 


In equation (3) m is a proportionality constant 
which will depend on the pecularities of the 
system used. The r—c relations obtained 
experimentally have been analysed with the aid 
of equation (3). Plotting r versus ¢ and taking 
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equal increments Ac, the corresponding r values 

can be used to find the magnitude of r_. The 

slope of the corresponding log (r — r_) — ¢ line 
m 

gives the value of aN’ 

It has been found that all the cases investigated 
can be represented by means of equation (3). 
Before giving the condensed results, as an example 
in Table 2 the r — c data for the filtration of Pur II 
over twillfill cloth are tabulated and in Fig. 4 
the r — c-line and the resulting log (r —r.) —e 
line are given. In Table 2: R, Ry_,o- Large 
values of R,, R, indicate that in these measure- 
ments the cloth has been clogged appreciably, 
giving too large values for r. The points indicated 
by (1) and (2) (see Fig. 4 and Table 2) have 
therefore not been used in drawing and using the 
r—c-line. In Table 3 the collected data are 
given. For the systems indicated by + no 
blocking filtration has been found. In these 
eases the value of (rz), ,o and of al could 
therefore not be evaluated. The value of 1, 
given here is the resistance of the cake at the 
lowest concentration c, used. In case (11) the 
blocking filtration has not been measured and 


1 1:9 2-0 
3 | 2-3 3-0 
4 1-9 2-0 
| 
5 20 | 1-7 
12 | 2-2 1-8 ; 
| 5 
13 | 40 hs 
2-1 | 16 
VOL, 
6 
956/57 
18 (2) | 1-2 70 
| | 
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Table 3 


| | 


Substance Filter | 


Pur I | 14000 
Pur I 


+ Glass p. cale 


aN 


Purl A 3640 15 
Pur 1B 140 +10 
Pur 28 
Pur 1D | cale : 

710-5 
Pur ILA 3:3 
Pur ITB . 5 6-6 

+ Pur Il 


Alum. 
Alum. 


103 
10°3 


aN 


| 


| 18680 36 10% 

5600 33 104 
810 3-104 

10°53 | 


3360 18 
1260 | 70 108 


T700 49 10% 
2-104 


105 
105 
105 
105 
105 


0-9 105 
10 105 
1-1 108 


° 
log (r-rao) 


Fic. 4. The r —ec and log(r —r,) relation for 


Puritine II filtered over twillfill cloth. 


c, and r, here have the same significance as 
mentioned just before. 

For the glass pearls the value of aN could be 
calculated from the size analysis of the product 
[1]. In the series (1) and (2) uncertainty existed 
over the value of r, as measured. These values 


have not been included in the evaluation of 


m 
aN 
aid of 


The r, as given has been calculated with the 


m 
v found from the other points measured 


for the determined value of c, given. The value 


of m — has been calculated from = and aN 


S N Ss 
by multiplication. 

Finally some additional informations on Norit 
and its three fractions are collected in 
Table 4, in which the values of r, determined 
at a filtration pressure of 0-1 atm and filtered 
over twillfill cloth are given. 


sieve 


Table 4 


| >= 10-10°5/ ¢ = 1-10°5 
Norit 871 
Norit A 222 
Norit B 417 
Norit C 838 
83-5 percent A 


+ 16-5 percent 
(weight basis) | 
+ 8 per cent C | 


(weight basis) | 


566 
92 per cent A 
420 


be drawn. 
The relation between r and ¢ is given by equation 
(3). The values of aN are such that at a not 


The following conclusions can 
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| 0) 
16 
(4) 
iv? 2-7 
(9) 210+ ‘ 
(6) ‘ 
7-10 
—5 
1-10 3-4 
(8) 104 | 
2°: 0 
(9) 66 
( | 
(10) 110-5 | 825 | 180 | 
0-40 1¢ 770 185 18 
1-10 | 560 | 185 1-3 
1-10°5 535 | 17 
70 | 10 on 
(11) Pur Il I 0-10 calc : 1-10-5 | 585 | 10 
96 10 ) 
0-16 80 570 | 280 72 195¢ 
| 83 10% 75 |2 
(12) I 0 | 435 | 
10 135 | +10 | 190 
(13) I » 3.10-3 
0-10 | 200 | 310 7-104 
4 
4 
4 
= 4 = 
| 
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very large concentration the value of r differs 
less than the experimental error from r_. As 
an example may be given that if the experimental 
error in r is assumed to be 1 per cent, for series 
(8) this concentration is at 10-5 10°, for series 
(7) at 8-8 10°. Therefore it is comprehensible 
that r rapidly reaches an apparently constant 
value ; r, can be taken for this value. 

The value of r, ~ (r4)g,9 a8 found earlier for 
one system. This fact is not surprising. For 


the systems investigated rn is constant (see later). 


Therefore (74), ,9 is independent of the concentra- 
tion. Increasing the concentration shows at 
c, the gradual transition from blocking filtration 
to cake filtration. No abrupt change of the 
corresponding resistances can be expected of 
the very first layers of solid substance deposited. 
In agreement with paragraph (a), the value of 
r. is independent of the filtration pressure, 
whereas at a lower concentration an influence of 
the pressure manifests itself in the influence of 
the initial filtration velocity, which governs the 
packing of the cake. 

The parameter m introduced can not be 
calculated directly because — is unknown. 


However in reasonable approximation for the 
cases investigated may be taken equal to 


unity,* in which case m can be read directly 


*In the publication on blocking filtration [2] it has 
been said that although in general 0 + S, in the cases of 
substances having a rather wide size spectrum containing 
a large amount of relatively small particles, O was equal 
to S. 

This conclusion was based on the fact that for these 


Oa 
substances fi r different concentrations - was found to 


be constant. The reasoning behind this was, that O and S 
change with the same variables but that in general O 
is also influenced by the form of the pores and of the 
particles. If however for one system the influence of 


OaNn 


the concentration is the same for S and O is constant 


and therefore also - it seems likely that the extra influence 


of the form of the pores and the particles is not very large 
and therefore 0 ~ S in these cases. However for the sieve 


fractions investigated (Exp. 4, 5 and 8), the assumption 


from Table 3. It can be concluded that m shows 
a great difference for Pur I and its fractions for 
which m ~ 2 and for the coarser Pur II and its 
fractions for which m ~ 75 filtered over nylon. 
It may be remarked that the values found for 
(11) may be too high, because they have been 
calculated with the aid of r, as given which 
may be lower than the resistance at the transition 
point which has not been determined. The same 
holds for (7). A very great influence of the 
particle size spectrum and of the absolute size 
of the particles on the value of r, is found. 
The influence on r_ is far less, if present at all. 

The impression is gained that in the close 
packed cakes especially the smallest particles 
give a great increase of the resistance ; compare 
(1) with (7) and with (4) and (5). This is not 
surprising. The denser a bed is packed, the 
larger will be the clogging effect by the smaller 
particles of the voids between the larger particles. 

The strong influence of the small particles is 
again confirmed by the experiments with Norit 
and its three sieve fractions. A more systematic 
research on the influence of the size-spectrum 
on the cake resistance as well as on the significance 
of the parameter m is evidently necessary. 


that 5 1, seems far less justified, even if it is keot in 
mind that these fractions are far from monodispers. It 
had namely been found for the sieve fractions treated in the 


OaNn 
sume paper that =o «(was not constant at different 


concentrations. However as has also been observed the 
values of S found for the lowest concentration used were 
always smaller than 1, which is a physical impossibility. 
All the values of S for these fractions can now be corrected 
by adding up a constant amount to each, such that not 
OaNn 
only ; becomes constant for each concentration, 
but also that each S is larger than 1. The calculated 


values of aN - taken equal to 1], are given below. 


The order of magnitude of these values is the same as 


those found by calculation with a mean value of a based 


on an apparent density of puritine of 1-0 and N 2250. 
60 — Wy: 11-10, (60 — 90) + (90 — 1504p: 
90 150 p: 21-10°5, (60 — 90 


(150 220): 37-105, 150 — 2204: 27-10°, 


. 
. 
i= 
q 
| 
+ 
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CONCLUSION 


For the systems investigated and for rather 
thin cakes it has been found that the resistance 
of the cake is determined, besides by the nature 
of the system, by the initial filtration velocity 
and by the concentration of the suspension. 
Both, by increasing give rise to a decrease of 
the specific resistance. For a constant initial 
velocity the relation between r and c can be given 
by the equation : 

mt 

in which r, is the resistance at the concentration 
c, and r_ is the asymptotic value of the resistance 
at an infinite concentration. The size and size 
spectrum have a large influence on the magnitude 
of r,, The small particles determine the magnitude 
of r, to a high degree. Many and relatively small 
particles increase the value of r,. 


[1] Heertes P. M. and van per Haas H. Studies in filtration, 
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A graphical method for the calculation of cooler-condensers 


G. H. P. Bras 


Chemical Construction Corporation, New York, U.S.A. 


( Received October 1956) 


Abstract—The method of Cotsunn anp Hovucen for the design of cooler-condensers is 
simplified by the use of a new graphical procedure for the determination of the liquid-gas interface 
conditions, by which trial and error computations are entirely eliminated. 


The design method is based on the application of a diagram of the logarithm of the non- 
condensing gas partial pressure versus the temperature. 


An example is given to illustrate the design method. 


Résumé— L'auteur présente une simplification de la méthode de CoLpuRN et HouGEN pour le 
caleul des condenseurs réfrigérants, par utilisation d'une méthode graphique pour la détérmination 
des conditicns & l'interface éliminant ainsi tout tatonnement et cause d’erreur. 


La méthode graphique est basée sur représentation du logarithme de la pression particlle du 
gaz non condensable du mélange gaz-vapeur saturé en fonction de la température. 


Tue method of and Hovucen [1] for 
the design of cooler-condensers for saturated 
gas-vapor mixtures requires the point to point 
solution of the equation, 


he (t; t.) h (t, t,) 4 k, Ai (Pp — Pi) 
U (t, — t,). (1) 


The heat transfer coctlicient, hy, is a combined 
coctlicient for the cooling water side. the tube 
wall, scale and dirt deposits, and the film of 
condensate. It includes all heat resistances 
except those at the gas side, and it can easily be 
calculated by conventional methods. For most 
engineering calculations the combined heat 
transfer coeflicient, ky, varies little in the cooler- 
condenser, and can be taken as constant. 

The gas-side heat and mass transfer coeflicients 
h, and k,, can be calculated by any of the available 
correlations in the literature, e.g... by 


and Co.Burn’s [2] equations, 
eG 

Pri 

JGM, 


3 
Pe M,, Sci (3) 


L’auteur donne un exemple pour illustrer cette méthode graphique. 


If the inlet and outlet temperatures for the 
gas-vapor mixture and the cooling water are 
fixed, the intermediate values for the cooling water 
temperature, ¢ can 


» and the gas temperature, t 


be found from a simple heat balance in nA of 
saturated gas-vapor mixtures, provided the gas- 
vapor mixture remains saturated during its flow 
through the cooler-condenser. 

Inspection of the equations (1), (2) and (3) 
reveals that equation (1) contains two unknowns, 
namely the temperature and partial vapor 
pressure at the gas-liquid interface, t; and p,, 
while equation (3) contains one unknown, 
namely p,;, because (p, p)/\n [(P — p,) 
(P — p, )]. 

Co_Burn and Hovucen [1] have therefore solved 
equation (1) by trial and error. After the calcula- 
tion of various intermediate values of the total 
heat flux, U (t, 


be calculated by graphical integration of the 


t.), the tube area required can 


equation, 


dQ = U (t, —t,) dA (4) 


q 
| 
(2) 
or, 
dQ 
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Equation (1) neglects to account for the sensible 
heat of the diffusing vapors. A correction factor 
to take care of the sensible heat of the diffusing 
vapors was first proposed by AcKERMANN [3]. 
Based on this work, the author [4] has later 
derived equations for the apparent heat transfer 
coefficients at the gas side, and at the condensate 
side of the gas film, respectively. By the applica- 
tion of the correction factors involved, equation 
(1) can be more generally written, 


ho (t; — 
ee 
h,. (t, —t,) +k, Ai(p, —p,) (6) 
or, 
hy (t; — 
. (t + hy — py). (7) 


The numerical value of the correction factors 
in equations (6) and (7) is usually near to unity 
in most industrial cooler condensers. Thus, the 
correction factors may be neglected, except under 
extreme conditions. For the present purpose, it 
will be assumed that equation (1) is sufficiently 
accurate. 

The trial and error calculations required for the 
solution of equation (1) are very time consuming. 
The author [5] has therefore proposed to solve 
this equation by a graphical procedure in a 
temperature-vapor pressure diagram. 

Cairns [6] has proposed to express the relation 
between heat flux, U At, and heat transferred, 
Q, by the equation for a parabola. This enables 
the calculation of the tube surface by analytical 
integration, if the heat flux is known for three 
values of Q. Cairns has obtained reliable results 
in four out of five cases. The author [7] has 
recently given another cooler condenser example, 
and the application of this design method also 
resulted in a relatively large deviation. 

In the same article [7] a further simplification 
to the Co_surn-HovuGcen design method was 
proposed. It was shown that the introduction of 
the group F’ = (L,/C,,). (Pr/Sc)*/*, permits the 
elimination of the mass transfer coefficient, ky. 
The factor F’ has the dimension of a temperature 
differential. It varies but slightly for most gas- 
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vapor mixtures over rather wide temperature 
ranges, so that it can often be taken as constant 
for design purposes. 

In a later development of the graphical design 
method for cooler-condensers, the author [8] 
has shown that the graphical procedure can often 
be shortened by the application of a polar diagram. 
The polar diagram cannot be of much help when 
the combined heat transfer coeflicient, hy, shows 
great variations between inlet and outlet of the 
cooler-condenser, however. 

There is thus room for further improvement in 
the design methods for cooler-condensers proposed 
so far. It is the purpose of this paper to present 
a new approach to the graphical determination 
of the gas-liquid interface conditions. This 
approach will eliminate all trial and error calcula- 
tions and will deliver the interface conditions 
within the accuracy of the graphical procedure in 


a very simple way. 


EQUATIONS 


From equation 3 follows, 


» 
k, ke -In (F Ps) (8) 
(Pre Pa (P Pr) 
where, 
iGs 
k je (9) 


M,, Sc 
If equation (8) is substituted in equation (1), 
the following relation is obtained after rearranging, 


hh (t; — ty) = hy (ty — ty) 
+ k,’ A{In(P — p,) — In(P — p,)} (10) 


or, 
— +h 4 
k kA 
in(P — p,) —p) (11) 
h, h, 
From equations (2) and (9) follows, 
= F. (12) 
h, ( m Se 


Substitution of R = hgh, and F =k, Aj/h, in 
equation (11) delivers, 


(R +1)t,— Rt, —t, + Filn(P — p,) 
F in (P — p,) (13) 
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A graphical method for the 


Next, the temperature t, will be defined, 


(R + 1) 
Substitution of equation (14) in equation (18) 
gives after rearrangement, 
R+1 
( F (t; t4) 


(14) 


In (P — p,). (15) 
When no use is made of the simplification 
obtained by the introduction of the Factor 
F =k,' A,/h,, equation (15) will appear in the 
form, 


te 
or, 
ho +h 
t4) = In(P — pj). (17) 


The equations (15-17) still 
unknowns ¢,; and p,. 


have the two 
However, the form of these 
equations enables the direct graphical determina- 
tion of the gas-liquid interface temperature, 
t,, by the following graphical procedure. 


GRAPHICAL PROCEDURE 

In Fig. 1, a plot is shown of In(P — p,) 
against the temperature. First, this plot can be 
used for convenient computation of the value 
for In(P — p,) as required for the calculation 
of t4, by equation (14). 

Equation (15), can be brought in the following 
general form, 


=In(P — p). (18) 

Equation (18) represents a straight line in a 
diagram of In (P — p) vs. t. The equation can be 
solved for two arbitrarily chosen values of t, 
e.g. tp and t,. This delivers two values of In 
(P — p), namely In(P — p)p and In(P — p),, 
respectively. The points D and E can be plotted 
in Fig. 1. A straight line drawn through the 
points D and E in Fig. 1 will intersect the curve 
for In (P — p,) in point C. Point C will have the 
coordinates In (P — p) = In(P -- p,), andt = t;. 

It can thus be seen that a direct solution of the 
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Fic. 1. Graphical determination of the interface 
conditions. 


interface temperature, t,;, is made possible by the 
use of a diagram of this type. 

Evidently, this method will be accurate within 
the limits of the graphical computation. For 
engineering purposes it is usually accurate enough, 
if the temperature ¢; can be taken the 
diagram within 0-1 or 0-2°F. 

Usually, it will be more convenient to bring 
equation (18) in the form, 


(R + 1) 
2-308 F 


from 


. (t — t4) = log (P — p). (18a) 
This equation enables to carry out the graphical 
procedure on commercially available semi- 
logarithmic graph paper, on which (P — p,) is 
plotted against the temperature. 


EXAMPLE 

In order to check the accuracy of this graphical 
procedure, the method has been applied to the 
design of a cooler-condenser for a saturated 
mixture of carbon dioxide and water vapor. 
Complete design calculations for this condenser 
have been published elsewhere [7], and reference 
is made to this publication for all details. 

The cooler condenser handles 7750 lb/hr of dry 
carbon dioxide, saturated with water vapor at 
latm. and 203°F and cools it to 95°F. The 
cooling water enters at 77°F and leaves at 149°F, 
and flows in countercurrent with the gas stream. 
The cooler condenser is of the shell and tube type 
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Table 1. Calculation of the reference temperature t 4 


| 

Point ty R R+1 ta 

No. °F *F °F Fr | °F 
1 203-0 149-0 2-84 2,506 0-001531 167-6 
2 202-1 141-75 2-98 2,510 0-001586 | 217-2 
3 201-2 135-6 3-12 2,515 0-001642 . 260-0 
4 199-4 126-05 3-39 2,516 0001746 — 322-9 
5 197-6 118-9 3-60 2,517 0-001830 373-5 
6 194-0 108-9 401 2,538 0-001976 443-7 
7 190-4 102-25 4-375 2,553 0-002106 ~ 487-7 
8 185-0 95-7 4-86 2,569 0-002280 528-5 
9 167-0 85-35 6-04 2,619 | oozes? | 581-4 
10 149-0 81-15 6-77 2,669 0-002910 604-0 
11 131-0 79-0 7-24 2,718 0-003030 618-7 
12 113-0 17-15 7-59 2,767 9003100 | 628-3 
13 95-0 j 770 | 7-00 2,803 0-003175 — 628-8 


and has 169 stainless steel tubes 1-05 inch O.D. points in the condenser and graphical integration 
and 0-824 inch I.D. on a triangular pitch. The of equation (5), where U (t, — t,,) = ho (t; — t,). 
gas flows at the shell side of the cooler condenser, First, the cooling water temperature, ¢,, is 
and the cooling water flows inside of the tubes. calculated for various intermediate temperatures 
The heat transfer coefficient, hy, is equal to of the gas-vapor mixture, ¢,, by the application 

71-8 B.t.u. /hr sq. ft °F. of a simple heat balance. If adiabatic conditions 

The calculation of the tube surface required in the cooler-condenser are assumed, and the heat 


% for the condenser, requires the computation of of the condensate is neglected, the following 
the heat flux, Ay (4; — ¢,,), at several intermediate equation applies for counter-current flow, 
: Table 2. Interface temperature and heat flux 

4 


Heat flux 


1 Via Graphical By Trial 
; | Procedure and error Deviation 
Point t, (t; — t,) ho (t; — ty) U (tg — be) 
No. °F F B.t.u./hr sq. ft. B.t.u./hr sq. ft. | - 
1 202-4 53-4 3830 3830 0-0 
2 201-4 59-65 290 | 4270 + 0-47 
3 200-2 64-6 4645 404) + O11 
+ 198-0 71-95 5160 5150 + 0-19 
5 195-7 76-8 5510 5510 0-0 
6 191-2 82-3 5905 5895 + O17 
7 186-3 84-05 6040 6025 + O25 
8 179-0 83-3 5980 5975 + 0-08 
9 153-8 68°45 4910 4925 0-30 
10 130-5 49°35 3540 3540 00 
1! 110°3 31-3 2247 2275 
12 95-5 17-75 1275 1263 + 1-00 
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Go di, L dt 
enthalpy of the saturated carbon dioxide-water 


In this equation, i, is the 


vapor mixture, B.t.u./lb of vapour free gas. 
Next, the heat transfer coeflicient, hy is 
calculated from point to point by equation (2). 
Values for R = hg/h,gand F’ = (L, C,,). (Pr Sc)! 
are also calculated from point to point. Inter- 
mediate values for the reference temperature, 
then be 


can calculated by equation (14). 


In this example, the term (P — p,) in equation 
(14) has been multiplied by 10, in order to avoid 
negative numbers for the logarithms of small 
numbers. Accordingly, the graphical determina- 
tion of the interface temperature, t,, is carried 
out in a diagram of In [10 (P 


temperature, see Fig. 2. 


p)| vs. the 


The course of the calculations for the first 


point (gas entrance) is now as follows : 
Gas temperature, t, = 203-0°F, 
Cooling water temperature, ¢ 149-0°F. 


0-8341 


to 
Partial water vapour pressure, p, - 
a.t.m. 

0-8341 18 
0-1659 44 
3-056 


Humidity, 2-056 |b /Ib. 


23,650lb_ hr. 
Maximum free area for gas flow at the baffled 


Total gas flow, 7750 


shell side, 1-755 sq. ft. 


Gas Mass velocity, 13,500 lb. hr 


*T55 
sq. ft. 
Viscosity of gas mixture, « = 0-0334 lb/hr ft. 
Heat capacity of gas mixture, ¢ = 0-392 


B.t.u. /Ib °F. 
Thermal conductivity, k 
ft °F. 
Prandtl number, Pr 
Prt = 1-047. 
Reynolds number Re 
105 18,500 
12 x 0-0334 — 
0-0050 (Reference 


0-01225 B.t.u./hr 


1-070. 


Heat transfer factor, j 
[2], Fig. 2). 
Heat transfer coeflicient, 
0-0050 0-392 x 138,500 
047 
- 25-30 B.t.u., hr sq. ft. “F. 
R = ho/h, = 71-8 /25-30 = 2-84 


h 
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F’ = (L,/C,,) 
reference [7}). 

In [10 (P — p,)] = In 1-659 

Reference temperature, 


(Pr Sc)3 2506°F (See 


= 0-5062 


, (2-84 x 149) + 203 — (2506 x 0-5062) 
167-6°F. 
Equation (18) gives, for tp, 170°F, 


3-84 
In [10 (P — p)] (170 + 167-6) 

2506 
for ley 


3-84 


0-516; 220°F, 


In [10(P — p)] . (220 + 167-6) 
2506 


0-593. 


The two points D, and E, are plotted in 
Fig. 2, 
D, and E, intersects the saturated curve 
in point C,, 


and a straight line drawn through 
representing the interface 


conditions. The interface temperature is 
read from the graph, or t; 202-4°F. 
Heat flux, ho (t; — t,,.) = 71-8 (202-4 — 149-0) 

3830 B.t.u./hr sq. ft.°F. 
The calculations are repeated for the other 


intermediate points in the cooler-condenser. 


In table 1, point to point values for t,, R, F’, 
(R +1) and are published 
previously [7]. The values for t, shown in this 
table have been calculated as shown above. 


In Fig. 2, 


given as 


the graphical determination of the 
interface conditions is shown (points C, — C\,). 
The abscissa of these points represent the interface 


Saturation curve,+ 
(P-p.). 10, vs.f 


Ln 
+ 


150 200 250 
Temperoture 


Pig. 2. Graphical solution of cooler condenser example. 
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temperature t;. The values for t; found graphically 
are listed in Table 2. 

Table 2 also gives values for (t; 
heat flux, hy (¢; —¢,). Further, point to point 
values of the heat flux calculated by trial and 
error are shown in this table as reported pre- 
In the last column, the percent 


and the 


viously [7]. 


deviation between the two values is shown. 


The deviations are due to the graphical computa- 
tion and they are well within the accuracy usually 


required for engineering calculations of this type. 


NOTATIONS 


A = Tube surface area, sq. ft. 


C. Cp» = Specific heat at constant pressure of gas-vapor 
mixture and vapor, respectively, B.t.u./Ib °F. 

Specific heat of cooling liquid, B.t.u./Ib’F. 
Cy, = Specific heat at constant pressure of gas vapor 


mixture, B.t.u./Ib mole °F. 

D = Diffusion coefficient, sq. ft/hr. 

= ky Ay/hy, °F 

F (L, Cc.) °F. 

G = Mass velocity of gas vapor mixture, lb/hr 
sq. ft. 

Go = Flow rate of vapour-free gas, Ib. /hr. 

hy = Combined heat transfer coefficient for cooling 


" water, tube wall, scale, and condensate, 

q B.t.u./hr sq. ft °F. 

a h, = Gas film heat transfer coefficient, B.t.u./hr 

3 sq. ft °F. 

Ai = Enthalpy difference between the vapor in the 
gas and the condensate on the tubes, B.t.u./Ib. 

of vapor. 
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j = Heat and mass transfer factor, dimensionless. 
k = Thermal conductivity of gas-vapor mixture, 
B.t.u./hr ft °F. 
k, = Gas film mass transfer coefficient, lb/hr. sq. 
ft. atm. 
Ky Peg, lb/hr sq. ft. 
L = Flow rate of cooling liquid, Ib. /hr. 


L,;, Lp = Latent heat of condensation, respectively, at 
t;, and at ¢,, B.t.u./Ib mole. 
M.,, M,, = Molecular weight of vapor and of gas-vapor 


mixture, respectively, Ib/Ib mole. 
P = Total pressure, atm. 
Pr Prandtl number, dimensionless. 
p, = Partial vapor pressure at saturation, atm. 
pressure in the 
the gas-liquid 


gas-vapor 
interface, 


Po Pi Partial 
mixture, and at 
respectively, atm. 

Ap = Pp — Py atm. 
Pgs = Logarithmic mean partial pressure of non- 
condensing gas, atm. 
Q = Total heat transferred, B.t.u./hr. 
R = ho/hg, 
Se = Schmidt number, y/p D, dimensionless. 
t = Temperature, °F. 
t4 = Reference temperature, see equation 14, °F. 


vapor 


dimensionless. 


tg, ti, Qe = Temperature of gas mixture, of gas-liquid 
interface, and of cooling water, respectively, 
F. 
At = t, —t,, °F. 
U = Overall heat transfer coefficient, B.t.u./hr 
sq. ft °F. 
‘= Ap h,. 
« = ky Ap dimensionless. 
Ay Ap = Latent heat of condensation, respectively at 


t;, and at ¢,, B.t.u./Ib. 
p = Density of gas-vapor mixture, Ib/cu. ft. 
= Viscosity of gas-vapor mixture, lb/hr ft. 
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Heat-transfer in cylindrical channels involving potential velocity 
distributions 
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Exact solution are presented for the problem of heat-transfer involving potential 


velocity distribution, in cylindrical channels, without making the usual assumption that 


1 oT 
r or 


ra As* r 

parameters |-,-,——, -| and not on | -, 
8 a ss 


the problem. 


Résumé 


Les auteurs preséntent des solutions exactes du transfert de chaleur avec distribitution 


The temperature distribution is found to depend on the dimensionless 


y*3 as in the approximate treatment of 
«@ 


de vitesses, dans les canaux cylindriques sans faire Vhypothése habituellé : 


YT 1 dT 
4 am 


* 


r 2 Vma As* 


r a 
et non de — 


8 8 a x 


l1 INTRODUCTION 


Torrer [1] investigated the temperature distri- 
bution for plug flow, in cylindrical channels, 
kept at constant wall temperature, with uniform 
and constant generation of heat within the 
fluid, when the fluid enters the channel at a 
prescribed temperature. He also considered the 
problem of temperature equalization, when the 
inlet temperature has one constant value for 
part of the radius, and another value for the 
rest of the radius. The solution of the first 
problem is helpful in estimating temperatures 
developed in chemical and nuclear reactors, 
and that of the second problem in estimating 
average life of non-uniformities of the temperature 
in such streams. 


His treatment of these problems is approximate 
and is based on the assumption that 
YT 
* yw 


The solution of the 


Ils trouvant que la distribution de la température dépend des paramétres sons dimensions : 


Comme dans le calcul approximatif du probleme. 


2838 


r ow 


x As* 


s sVm s «4 


simplified energy equation reveals that the tem- 
perature distribution depends only on _ non- 
dimensional quantities (r/s, 4s*/«) and a dimen- 


a 
sionless group 
s m * 
In this paper, we give exact solution of the 
energy equation corresponding to ToppeEr’s 
problem. We find that the temperature distribu- 


tion depends on the dimensionless quantities 
r « Viz 


, |. The dimensionless group 


8 ax 


a 
(; yo ") does not appear in the exact solution. 


2. Untrorm Heat GENERATION 


Consider a stream of fluid entering the channel 
of radius r = s. At 2 = 0, the stream tempera- 
ture is 7), and the uniform wall temperature is 
T,. The energy equation for steady state, and 
in case of axial symmetry simplifies as 


e 
‘ 
4 
Pa : 
: 
| | 
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@ is defined by the relation 


with the boundary conditions 


In case of potential flow, V. 
and (1) becomes 


The complete solution of equation (5) is 


nel 
is the Bessel 


where B, and 8 are constants. 


The condition (3) is satisfied if 8 


(8,), and (4) gives 


\’ J, ( w) 


Multiplying (8) by 


integrating within the limits zero to 1, we have 


3. TEMPERATURE MIxING BETWEEN 
CONCENTRIC STREAMS 


In this case, the temperature of the entering 
fluid (at 2 = 0) is assumed to be 7, for part 
of the radius, say from r = 0 to r =a, and 7, 
for the rest of the radius (from r =a to r = s). 
The wall temperature remains constant T,. 
The energy equation (1) with = 0 simplifies 
to 


row 

where T = 6 + T,. 

The boundary conditions are 

atz =0,0=0,=T,-—T,a>r>0, (12) 
a, 

at r=ae, = 0. (13) 


The solution is 


0 dD) Bae vans) 
nel 


are the positive roots of J, (8,) 0. 


(8,2): (14) 


where £8, 
B. are determined from (12). 


B, J, (8, w) J, (8,, w) w . dw 


0, | +0, | J, (8, w) wdw. (15) 
a oO m b ‘ m 


. 
a 


Hence B, 


B, iJ, (B,)}" 
a a 


a +- +. 
T + T., (2) 

@ = 0, for z > 0 and r = s. (3) 

= 0, = T, — T,, for x = 0 and r < s. (4) 

a 

190, ,, A e* 

+ —¥0 (5) ; 
where A . 2x). (6) 
VOl 
6 
195¢ 
4 function of the first kind and of zero order. ae 
4 
| 
+ (1 — w*) (8) 
4 nel a/s i 
where w = (rs). 
(8,, 0) w dw, where £,, 
need not be the same eigen value as and 
20 As? 1 
q B, J; (8,,) | | 
; 
By? T, — 7, B, (B,)} | 
nel nel 
Vin (axe aT,—T, a 
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4. Discussion OF THE RESULTS 


From equations (10) and (17), it is obvious 
that the temperature of the fluid at any point 
is a function of 

ss a 
Toprer’s example, the temperature depends on 


whereas in 


(’ plots the temperature distri- 
ssV, 8s «@ 


bution on the axis of the channel (r = 0), against 
the abscissa (a7 +). Here it is not possible 
4 m 
to compare the results with those of Topper’s, 


because the temperature cannot be plotted as 


a function of (5 . 


4. NOTATIONS 


A = (q/pc), uniform heat generation. 
Jo, J, = Bessel function of the first kind, and 
of zero and first order. 


B,, = coefficient of characteristic function 
in various equations, 
temperature, initial temperature at 
initial temperature at 
initial temperature, wall 
temperature. 
mean velocity, local axial velocity. 
= radial coordinate at discontinuity of 
initial temperature. 
= heat capacity. 
heat generation per unit volume and 
time. 
radial coordinate. 
= tube radius. 
= (r/s). 
axial coordinate. 
thermal diffusivity. 
B, = @ zero of Jo. 
p = density. 
Gat w = 0, and = (T — T,)e™. 
A — (V,,,/22). 


r< a, 


In conclusion, I thank Dr. S. D. Nigam for 
suggesting this problem to me. 
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Vicror Goin: pH Measurements. Methuen Mono- 
graph, 1956. 125 pp. 9s. 6d. 


Tare is a story of a Professor who remarked, on readin? 
a student's doctoral thesis, * There is much that is true 
and much that is new in this thesis, but that which is 
true is not new, and worse, that which is new is not true.” 
Dr. Goin, in his Methuen Monograph on pH measure- 
ments, has tried hard to avoid the second fault, and has 
succeeded completely in keeping to established theoretical 
and experimental results. 


W here, however, there is less certainty about the measur- 
ing of pH, as for example near detergent micelles, polyelect- 
rolyte molecules or charged surfaces, Dr. Gown offers no 
suggestions, preferring to remain on the solid ground of 
long established fact rather than venture into the quick- 
sands of contemporary research. 


An author's approach, however, is always a matter of 
taste, and the book must also be judged on the material 
therein as well as on its general scope. There is no doubt 
that Dr. Goin has been meticulously painstaking in 
selecting his pH standards and in discussing how they 
should be used. The book is thus of interest to anyone 
working on very accurate pH measurements (i.e. to better 
than 0-01 pH unit) and to those interested in the physical 


meaning of pH determinations. 


In the reviewer's opinion, however, this is a book on pH 
which virtually ends the fundamental importance of the 
pH scale : now that there are very fundamental difficulties 
in reconciling the different definitions and measurements 
of pH (i.e. peH, pH, pwH, psH), might not anyone working 
outside the field of acid-base catalysis prefer to quote 
the result of a measurement with some standard pH 


meter and leave it at that ’ 


Where a fuller discussion might have been useful is in 
the field of non-< queous solvents. The recent importance of 
titrating long-chain acids and polymers in non-aqueous 


solvents would seem to have justified this 


Finally, who is the book intended for’ The research 
physical chemist would probably turn up the fourth 
edition of Britton’s book “ Hydrogen Ions,” (1956), 
or Bares’ book “ Electrometric pH Determinations ” 
(1952). The chemical engineer would, like the reviewer, 
probably be content to read the instructions on the pH 
meter. However, the industrial research chemist who 
requires a ready guide to pH determinations, who requires 
considerable accuracy or who wants to measure a pll 
which has some reasonable physical interpretation, will 
find this little book a clear, useful and accurate bench 


companion, 
J. T. Davies 
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Physical Methods in Chemical Analysis. Vol. III 
Edited by W. G. Bert, Academic Books Ltd., London, 
1956. xii + 652 pp. $15.00. 


To AN increasing extent the science of chemistry is being 
facilitated by the development and availability of physical 

methods for qualitative and quantitative analysis in the 
widest sense of the word. Books concerned with these 
newer specialized methods are becoming more prevalent ; 
in essence they must provide reasonably comprehensive, 
clearly-written reviews that are understandable to the 
non-expert. In general, Bert's Volume III maintains the 
high standard of the earlier volumes, and some of its 12 
chapters are excellent. In others the Reviewer feels that 
the space allotted has been too small so that, although the 
theory and general principles are stressed, and rightly so, 
the number of actual examples of the scope and uses of the 
particular technique is too small. A book of the same length 
covering 6-9 of the topics would probably be of more 
value to the chemist, chemical engineer, or analyst. 
Two of the topics discussed — gas chromatography and 
nuclear magnetic resonance spectroscopy are of par 

ticular current interest. It is only a few years since the 
potentialities of these techniques were realized, but today 

relatively inexpensive gas chromatography instruments 
and distinctly expensive, but none-the-less highly useful, 
nuclear magnetic resonance spectroscopy instruments are 

available commercially. The chapter on gas chromato- 
graphy, written by a pioneer of the technique, provides 
an excellent introduction to the subject, but is dis 

appointingly short (27 pp.) compared with many of the 
other chapters. Thus the reader must seek elsewhere 
for more detailed information, and, although the list of 
35 references includes most of importance up to 1955, 
one wishes that it were more comprehensive. The section 
on nuclear magnetic resonance spectroscopy 1s longer, 
with an adequate reference list. Written by a leader in 
the field, it constitutes one of the best reviews yet available 
and covers theory, experimental methods, and the inter 

pretation and uses of this important branch of spectros 
copy. The uses of neutron spectroscopy in chemical 
analysis are admirably summarized, often in tabular form, 
in the longest and most comprehensive chapter of the 
book. Electrochromatograply is another relatively 
new technique and its practical applications, particularly 
in the biochemical field, are emphasized in a review that 
sets forth the principles and describes in some detail the 
apparatus required; the diagrams in this chapter are 
party ularly good. Other topics covered are electroanaly- 
tical methods in trace analysis, the high frequency method 
of analysis, field emission spectrose »py, flame photo- 
metry, microwave spectroscopy, flucrescent X-ray spectro 
metric analysis, and analytical distillation ; a section on 
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the theory and principles of sampling for chemical analysis 
is a useful addition. 


The book is well printed and ressonably free from errors. 
Despite the high price that one has come to expect for a 
compendium of this type, this is a recommended book 
with which all practising chemists and chemical engineers 


should be familiar. . 
R. N. 


McCane and Surra: Unit Operations of Chemical 
Engineering. McGraw-Hill, New York, 1956, 956 pp. 79s. 


Tus book has been written for students who are meeting 
unit operations for the first time. It assumes only a 
knowledge of elementary mathematics and physics. In 
common with previous books on the same subject, a broad 
interpretation is taken of the scope of unit operations. 
Over a quarter of the book is devoted to fluid mechanics 
and heat transfer. An excellent introduction is given to 
fluid friction and boundary layer theory. The importance 
and applications of dimensional analysis are stressed 
throughout the book. 


reviews 


All of the unit operations normally covered in a 
University course are discussed. Multicomponent separa- 
tions, transient phenomena, and operations at extremes of 
temperature, pressure, or velocity have been deliberately 
omitted. Some discussion of these points would be 
desirable, particularly since these fields are becoming 
increasingly important. 


The chief criticism that can be offered is that the book 
is too long, and, as a consequence, too expensive. Books 
intended for undergraduate use should be carefully pruned 
of all material not absolutely essential if they are to be of 
the greatest value. This might be done by eliminating 
the long introduction (which is essentially an introduction 
to chemical engineering), and by reducing the large amount 
of space given to description of equipment. For example, 
twenty-three pages are used to describe pipe, pipe fittings, 
and valves. The descriptions of equipment are, in fact, 
excellent, but a separate book might be a better place for 


most of them. 


Professors McCase and Smiru have, however, produced 

a very readable book which deserves a place on the 
chemical engineering bookshelf. 

G. A. Ratcuirr 
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SOME CURRENT PAPERS IN OTHER JOURNALS 
Chemical Engineering Progress 


VOLUME 538, No. 4 (partial contents) APRIL 1957 
A. Cnarnes and W. W. Cooper: Linear programming—solving management problems 
A. C. Dresurrecp, Jr. : Mechanism of hot surface drying of fibrous sheets 
F. G. Kerry : Safety aspects of air separation plant cycles 
E. C, CLark (Edited by) : Discussion of air separation plant safety 
E. Karwat : Linde’s view on air separation plant safety 
W. E. Karz: Status of electric membrane demineralization 
J. A. A. and W. Sommonps : Mixing patterns in helical-flight mixers 


N. M. Levitz, E. J. Perxus, H. M. Katz, and A. A. Jonxe: Fluidized bed process 
produces UF, 


P. A. Haas, E. O. Nurmi, M. E, Wuatvey, and J. R. Encet : Midget hydroclones remove 
micron particles 


ERRATUM 


Hovucurton G., McLean A. M., and Rrrenre P. D. Compressibility, fugacity, and water solubility of carbon dioxide in 
the region 0-36 atm and 0-100°. Chem. Engng. Sci. 1957 6 132. 


The first line of equation (3) in the left-hand column of page 133 should read : 
P 


= 
'* RT 
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